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ABSTRACT 



u: 
^' 

Of 

Q ■ We explore new observationally-constrained sub-resolution models of galactic out- 

\ flows and investigate their impact on the circumgalactic medium (CGM) in the redshift 

. range z = 2 — 4. We perform cosmological hydrodynamic simulations, including star 

formation, chemical enrichment, and four cases of SNe-driven outflows: no wind (NW), 
an energy-driven constant velocity wind (CW), a radially varying wind (RVWa) where 
the outflow velocity has a positive correlation with galactocentric distance (r), and a 
^ ■ RVW with additional dependence on halo mass (RVWb). Overall, we find that the 

CSJ , outflows expel metal-enriched gas away from galaxies, significantly quench the star 

OO ■ formation, reduce the central galactic metallicity and enrich the CGM. At z = 2, the 

I radial profiles of gas properties around galaxy centers are most sensitive to the choice 

, of the wind model for halo masses in the range (10^ — 10^^)Mq. 

■ We infer that the RVWb model is similar to the NW case, except that it substan- 

tially enriches the CGM: the carbon metallicity {Zc) is 10 times higher in RVWb than 
. in NW at r ^ and the warm gas of 10^ — 10^ K and 5 < 100 is enriched at a 

level ^ 50 times higher than in NW. Moreover, we find that the impact of models CW 
and RVWa are similar, with the following notable differences. RVWa causes a greater 
suppression of star formation rate at z ^ 5, and has a higher fraction of low-density 
' (S < 10), warm-hot (lO'' — 10^ K) gas than in CW. Zc in the CW model increases 

H . monotonically with S a.t 6 ^ 10, while in RVWa Zc first decreases and then remains 

constant over the same density range. Outflows in CW produce a higher and earlier 
enrichment of some IGM phases than RVWa: the warm 10"* — 10^ K, 5 ^ 10 — 100 
phase is enriched 10 times more and at z = 4 the underdense S ^ 0.1 — 0.8 warm-hot 
10^ — 10^ K gas has already a median Zc ~ 0.002Zc,q in CW. By visual inspection, we 
note that the RVWa model shows bipolar outflows and galactic disks more pronounced 
than in all the other wind models. We present simple fitting formulae for the Zc as a 
function of gas density and of the galactocentric distance, also for the abundance of 
triply ionized carbon (CIV) as a function of the galactocentric distance. We predict 
observational diagnostics to distinguish between different outfiow scenarios: Zc of the 
CGM gas at r ~ (30 — 300) ft." ^ kpc comoving, and CIV fraction of the inner gas at 
r < (4 — 5)h^^ kpc comoving. 

Key words: Cosmology: theory - Methods: Numerical - Galaxies: Intergalactic 
Medium - Galaxies: formation 
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1 INTRODUCTION 

Energy feedback from star formation (SF) and pow- 
erful supernovae (SNe) explosions are believed to 
eject _gas fro m gala xies driving outflows called 



20051 ). 



wnd jil (e.g., Ijohnson fc AxfordI 197ll: Mathews fc Baker! 



197ll: iLarsonl 1 19741 : IVeilleux. Cecil fc Bland-Hawthornl 



SNe-driven outflows are 
servations at low redshifts £e 



detect ed in many ob- 
Lvnds fc Sandagd 



19631: iMcCarthv. van Breugel fc HeckmanI 



Bland fc TuUv 



198 

Lehnert fc Heckman' 



.Heckman. Armus fc Milev 
199dl: IStrickland et all l2000l . 



Rupke. V cillcux & Sanders l2005l). reaching 
high as 1000 k m/s (IDiamond-Stanic et al 



at high- 2: (e.g., Franx et al.l 1 19971: iPettini et al 



veloc ity as 
also 
20011: 



201 



Adelbereer et al.l l2003l : IShaplev et alll2003l : iNewman et"all 
2012b|)^ extending up to di stances of 60 — 130 kpc physical 
fe.g.. lLundgren et al.l[20li ). 

We use the term circumgalactic medium (CGM) to de- 
note the gas-phase structures, excluding star - forming gas, 
within 300 kpc (physical) of galaxies, following ISteidel et al.l 
l|2010ll , from whom we derive our model. Galactic winds 
are considered to be the primary mechanism by which 
metals are ejected out of star-forming regions in galaxies 
and deposit ed into the CGM and the intergalactic medium 



(IGM) (e.g..lLarson fc Dinersteinlll975l : [Aguirre et"all 



Scannapieco, Ferrara fc Madau' g ' Aracil et al. ' 



200 1| : 

200J; 



Fox et al. 2007; Bouchc ct al. 2007; Kawata . 



Fieri. Martel fc Greno: 
20101 : iGauthier fc Che: 



3[2Qi3; 



2004, 
Ranch 20071: 



IPinsonneault. Martel fc Pieril 
I). The IGM can also be metal- 
enriched by active galactic nuclei (AGN) feedback (e.g. , 
iKhalatvan et al.ll200^ : ICermain. Barai fc Martej |2009l ') and 
heated by blazars (e.g.. Puchwein et al.,.2012 ). 

SNe and starburst driven outflow is an important source 
of feedback in galaxy evolution, and constitute a key in- 
gredient of current galaxy for mation studies, both i n hy- 
drodynamic simulatio ns (e.g.. ISpringel fc He rnguist 20031. 
hereafter SH0 3. .Schave et alj|2010l) and in se mi-analytical 
models (e.g., iBertone. Stoehr fc Whit3 bOOSl ). The winds 
act by removing gas available to make stars, hence quench 
SF, and are argued to suppress the formation of low- 
mass galaxies, flattening the low ma ss end of the lumi- 
nosity function in simulations (e.g.. [ Theuns et al.l |2002| : 



iRasera fc Tevssiei|[200^ : IStinson et al.ll2007l 'l. which is closer 
to observations. Feedback is invoked t o reproduc e the 
realistic disk galaxies (e.g.. IWeil. Eke fc Efstathi oul Il998l: 
Sommer-Larsen. Gotz fc Portinaril 200(^!~ Governato et al.l 



20041 : [Robertson et all |2004| : lOkamoto et all bOOSl ). These 
outflows are also argued to affect the large column den- 
sity parts of the Lyman-Q absorption line forest (LLSs, 
subDLAs and DLAs) seen in the spectra of distant 
quasars, which traces the IGM matter dist r ibution (e.g.. 



McE>qnald et al.l l2005l: iKoUmeier et al.1 l2006l : iTescari et al] 
20111 : IViel. Schave fc Boothll2012l ). 



The detailed internal physics underlying the ori- 
gin and driving of galactic outflows is complex, oc- 
curing on scales proper of the multi phase structure 
of the interstellar medium (ISM) (e.g., iHeckmanI l2003l : 



^ Here, we use the terms wind and outflow synonymously, mean- 
ing continuous outward flow of gas from a galaxy, which might or 
might not escape depending on its velocity and galactic potential. 



IStringer et ahl |2012| ). The gas is likely accele rated either 



by thermal pressure (e.g.. Chevalier & Clogg 1985), radi 



ation pressure (e.g.. iMurrav. Quataer t^ fc Thompson! [20051: 



Sharma. Nath fc Shchekinov 20l'il : Chattopadhvav et al.l 
2012! : IZhang fc T hompson"20lJ~ ram pressure , cosmic rays 



(e.g., ISamui. S ubramanian fc Srianandl |2010!: lUhUg et al l 
2012!): or a com bination of them (e.g.. iNath fc Silk 20o'9r 



Everett. Schiller fc Zweibel .201Q : ISharma fc Nath 201^; 



Cosmological simulations do not resolve the scales at which 
these physical processes actually happen. Outflows are hence 
implemented in a phenomenological way in the simulations 
using sub-resolution prescriptions and their interplay with 
the larger-scale environments is investigated. 

Energy ejection by SNe and starburst in galaxies, of- 
ten in the form of powerful outflows, can be incorporated in 
a number of ways as a source of feedback in the numer- 
ical scheme of cosmological simulations. A thermal feed- 
back, where SNe energy is distributed as heating energy 
of ne ighboring gas, is wel l-known to be ineffective (e.g. . 



Katz 1992: 'Fricdh fc Bcnz"l995': 'Steinmet z fc Mulled Il995l : 
Katz, Weinberg & Hcrnquist 1996), because the dense star- 
forming gas has high cooling rate, and the injected ther- 
mal energy is radiated away quickly before it can signif- 
icantly impact the gas. Therefore depositing the SNe en- 
ergy in the kinetic form is a more popular implemen- 
tation in the literature, which h as been shown to have 
signiflcant feedback effec t s (e.g.. iNavarro fc White! [l993l: 
Mihos fc Hernguist! '1994'- 'Cen fc OstrikeiJ '20001; 'Kawata 



2001!: ISpringcl & Hernguist 2003; Dalla Vccchia fc Schw 



2008! : iDubois fc Tevssieri 120081 : lOppen hcimcr ct al. I l20m 



In our current work we adopt this kinetic feedback by giv- 
ing a velocity kick to the affected gas. 

Alternative forms of SNe feedback have also been im- 
plemented for example by considering that the affected 
gas undergoes adiabatic evolution (e.g., Mori et al] Il997l : 
iThacker fc Couchman!l2000l : |BrooES^[200i); by turning 
off radia tive cooling temporarily for part of the neighborin g 
gas (e.g.. lGovernato et al.!l2007! ; |Piontek fc Steinmet3l201ll ): 
or by distributin g SNe energy to hot and cold gas phases 
separately (e.g., Marri fc White! 12003! : IScannapieco et aP 
I2OO6! : iMurante et al.! I2OI0! ). Note that few studies have 
proposed prescriptions fo r efficient thermal feedback 
in SPH s imulatio ns (e.g., Kav. Thomas fc Theuns! 12003! : 
iDafla Vecchia fc Schave!l2012! ). 

Our kinetic feedback models ( t|2.2l H2.3|l are based 
on the so-called energy-driven wind scenario origi- 
nally proposed by SH03, and subs eq uently used by 
others (e.g.. [ Tornatore et al.l 12004!: iNagamine et al.l 



2007!: !Dafla Vecchia fc Schave! !2008!: !Tescari et al.! !2009! : 



Fabian et all I2OI0! : I Barnes et al.! I201J ). Here, a fraction 
of SNe energy provides the outflow kinetic energy, and 
the wind speed is independent of galaxy mass or SF rate. 
If we consider the underlying physics, which occurs on 
scales orders of magnitude below the scales resolved in 
cosmological simulations, such an outflow is likely to be 
driven by the thermal pressure of SNe, and might be (more 
physically) called a thermally -driven wind. 

Other wind models in cosmological simulations gen- 
erally involve a different formulation of outflow velocity 
and wind mass loading factor (described more in >i2.2|l in 
terms of galaxy properties (mass, vel ocity dispersion, SFR), 
sometimes suggested by observations. lOppenheimer fc Dave! 
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ll2006t) implemente d mom entum-driven wind, (also e.g., 
iTescari et al.l '2009'. '2 01ll'l. using an alytical prescriptions 
from Murray. Q u ataert fc ThompsonI (i2005, ') for momentum 
injection provided by radiatio n pressure of photo n s and SNe 
(or, radiatively-driven wind); IChoi fc Nagamin^ l|201ll ) de- 
veloped a mul ticomponent and var i able v elocity galactic 
outflow model; IPuchwein fc Springell ('2012) implemented a 
halo mass dependent energy-driven outflow model. 

In this paper we explore new wind models, where the 
galactic outflow speed is a function of radius (distance 
from galaxy center ), as motivated by recent observations 
l|Steidel et al.l I2OI0I ). We investigate both the cases of a 
mass-independent dependence on radius, and a halo mass 
dependent parametrization of the radially varying wind. To 
our knowledge, this is the first numerical implementation of 
both of these outflow models. 

This paper is organised as follows: we describe our nu- 
merical code and simulation setup in 321 in 33] we present 
and discuss our results, while in i|4] we give a summary of 
the main findings and discuss possible future applications. 



2 NUMERICAL METHOD 

We use a modified version of the TreePM (parti- 
cle mesh) - SPH (smoothed particle hydrodynamics) 
code GADGET-3, whose first version was described in 
ISpringel. Yoshida fc W hTte (200t]. The public release of the 
code (GADGET-2, [Springcl 2005) contains a time inte- 
gration using energy and entropy c onserving formulation 
of SPH ISpringel fc Hernguistl [20o3 ). uses fully adaptive 
smoothing lengt hs, and a stand ard SPH artificial viscos- 
ity prescription (|Monagha3[l993). GADGET-3 includes a 
more efficient domain decomposition to improve the work- 
load balance over GADGET-2. 

Some of the additional subgrid physic^ included in the 
semi-public version of GADGET-3 code we use are out- 
lined in H2.ll The general galactic wind feedback model is 
described in i|2.2l our new outfiow models in H2.31 and the 
code implementations in i ]2.4l The simulations we perform 
are presented in ii2.5l 



2.1 Sub-resolution Physics: Radiative Processes, 
SF, Chemical Evolution 

Radiative coolin g and heating is based on the origina l im- 
plementation of iKatz. Weinberg fc Hernguistl (| 19961 ) and 
improved by adding metal-line cooling which is imple- 
mented by adopting t he cooling rates from the tables of 
LWic rsma, Schayc & S mithl l|2009h . Net cooling rates are 
computed element-by-element tracking 11 species: H, He, 
C, Ca, O, N, Ne, Mg, S, Si, Fe. A spatially-uniform time- 
dependent photoionizing background radiation is consid- 
ered &wn_the_cosmic microwave background (CMB) and 
the iHaardt fc Madaul (,2001) model for the ultraviolet /X- 
ray background produced by quasars and galax;ies. Contri- 
butions from the 11 elements are interpolated as a function 



^ By subgrid we mean sub-resolution, referring to physical pro- 
cesses occuring at length scales smaller than the resolved scales 
in our simulations. 



of density, temperature and redshift from tables that have 
been pre-computed using the public photoionizat ion code 
CLOUDY (last described bv iFerland et all 1 19981 '). assum- 
ing the gas to be dust free, optically thin and in (photo- 
)ionization equilibrium. 

SF and SNe feedback are implemented following the 
effective subresolution model by SH03. In this model, the 
physics of multiphase structure of the ISM, on scales un- 
resolved in cosmological simulations, is modeled using spa- 
tially averaged properties describing the medium on scales 
that are resolved. Gas particles with density above a limit- 
ing threshold, psF = 0.13 cm~^ (in units of number density 
of hydrogen atoms), are considered to contain two phases: 
cold condensed clouds that are in pressure equilibrium with 
ambient hot gas. Each gas particle represents a region of the 
ISM, where the cold clouds supply the material available for 
SF. Star particles are collisionless, and are spawned from gas 
particles unde rgoing SF, according to the stochasti c scheme 
introduced by iKatz. Weinberg fc Hernguistl (|l996l ) . We al- 
low a gas particle to spawn up to four generations of star 
particles; therefore a typical star particle mass is about one- 
fourth of the initial mass of gas particles. 

Stellar evolution and chemical enrichment feedback 
ar e implemented followin g the chemical evolution model 
of iTornatore et al.l |20o3). Production of 9 different metal 
species (C, Ca, O, N, Ne, Mg, S, Si, Fe) are accounted 
for using detailed yields from Type la SN (SN-Ia), Type 
II SN (SN-II), along with low and intermediate mass 
stars (LIMS) in the thermally pulsating asymptotic gi- 
ant branch (TP-AGB) phase. Contributions from both SN- 
Ia and SN-II to thermal feedback are considered. Mass- 
dependent time delays with which different stellar popu- 
lations relea se metals are included , adop ting the lifetime 
function by iPadovani fc Matteuccil (|l993l). Different stel- 
lar yie lds are used: for SN-Ia taken from Thielemann et ah! 
(l2003h. SN-II from IWooslev fc Wea^ jl995l) . and LIMS 



van den Hoek fc GroenewegenI ( 1997 ). The mass range 



from 

for SN-II is considered to be M/Mq > 8, while that for SN- 
Ia originating from binary systems is 0.8 < M/Mq < 8 with 
a binary fraction of 10%. 

We include a fixed stellar initi al mass function (IMF) 
according to the formalism given bv lChabriej l|2003l ). which 
is a power-law at M/Mq > 1 and has a log-normal form 
at masses below. However, we use power-law IMFs with dif- 
ferent slopes over the whole mass ran ge of 0.1 t o IOOM0, 
which mimics the log-normal form of Chabrieil (|2003l ) at 
lower masses, as tests indicate. In our model the functional 
form: 4>{M) — KM~^, is composed of 3 slopes and normal- 
izations: y = 0.2 and K = 0.497 for stellar masses 0.1 ^ 
M/Mq < 0.3, y = 0.8 and K = 0.241 for 0.3 < M/Mq < 1, 
and y = 1.3 and K = 0.241 for 1 < M/Mq < 100. Stars 
within a mass interval [0.1 — 40] Mq become SNe first before 
turning into black holes (BHs) at the end of their lives, while 
stars of mass > 40Mq are allowed to directly end in BHs. 

The chemical evolution model also incorporates mass 
loss through stellar winds and SNe explosions, which are self- 
consistently computed for a given IMF and lifetime function. 
A fraction of a star particle's mass is restored as diffuse gas 
during its evolution, and distributed to the surrounding gas 
particles. There is no AGN feedback in our simulations. 
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2.2 Galactic Wind Feedback 

We use subgrid models for kinetic feedback via SNe-driven 
galactic outflows, following the energy- driven prescription 
(© (originally from SH03), where the mass and energy car- 
ried away by outflows are regulated by two equations. The 
first equation relates the wind mass-loss rate Mw to the SF 
rate M*: 



(1) 



Here 77 is the wind mass loading factor, or the wind efficiency. 
Observations have shown that in galaxies outflow mass-loss 
rates are comparable to or few times l arger than SF rates 
(e.g.. lMartinlll999l: iPettini et al.1 12OO2I : iBouche et aH |2012| : 
iNewman et al. 2012a ). Thus, following SH03, we adopt a 
constant ?7 = 2 as the efficiency factor. 

The other equation relates the wind kinetic energy to a 
fixed fraction x of SNe energy; 

]^M^vl = x<^snM*. (2) 

Here is the wind velocity, esN is the average energy re- 
leased by SNe for each Mq of stars formed under the in- 
stantaneous recycling approximation. Combining Eqs. ([1]) 
and ((U, can be re- written as: 



1/2 



(3) 



Unlike SH03 . and following iTornatore et al.1 (j2007h and 
iTescari et all (|201ll '). we choose free parameter in 

our models, from which the effective x can be computed 
using Eq. Q. 

For our adopted Chabrier power-law IMF (^STTJ, esjv = 
1.1 X 10*" erg Mq^ . We initialize the wind energy fraction 
using = 400 km/s, which is the constant velocity of our 
simulation run CW ( i]2.5p . It corresponds to x = 0.29 of SNe 
energy being carried away by the wind. 

Outflow models implemented in the GADGET code 
generally involve different scaling relations of and 
rj in terms of galaxy velocity dispersion, mass, and/or 
SFR. For the original energy-driven outflows, Vm and 
ri are constant. The mome nt um-d ri ven wind formal- 
ism by lOppenheimer fc Pavel l|2006l . [ 20081) has oc 
(J ^JJLJIZ^H)--^ and rj cc 1/cr, where a is the galaxy ve- 
locity dispersion, and L/Lcrit is its lum i nosity in units 
of a critical value. In IChoi fc Nagamin3 ll201lD 's model, 
Vin OC SFR^^^ and ri is a f unction of galaxy stellar mass. 
IPuchwein fc Springell (I2OI2I ) assumed and r] to be pro- 
portional to halo mass. 

However, any dependence of on galactocentric ra- 
dius (or distance from galaxy center) has not been explored 
before and this is what we are going to do in the present 
work. 



2.3 Wind Velocity Dependent on Galactocentric 
Distance, v^ir) 

The new subgrid wind model is mo t ivated by the recent 
observational studies of ISteidel et al.l l|2010l ). Studying the 
metal-enriched gas kinematics in a region of ~ 125 kpc 
around star-forming (Lyman-break) galaxies at redshifts 
2 = 2 — 3, they are able to reproduce their spectroscopic 



data using a simple model for outflows and circumgalactic 
gas. Considering a spherically-symmetric gas flow, quanti- 
ties are parameterized as a function of galactocentric dis- 
tance, or radius r. The outflow acceleration is described as 
a power- law function of r: 



a{r) (X r " = 



dv{r) 
dt 



dv(r) 




— v{r) 


dv{r) 


dr 




dr 



(4) 



Such a p ower-law acceleration h ave also been adopted re- 
cently bv lCauthier fc Chen! l|2012l ) . Integrating this equation 
with appropriate boundary conditions one ca n get the out- 
ward v elocity v{r), which is given in Eq. (16) of lSteidel etHI 
l|2010l) . We adopt a similar expression for the radially- 
varying velocity of galactic winds, namely: 



(5) 



Here rmin is the distance from which the wind is launched 
and where the velocity is zero, ReB represents the outer 
edge of gas distribution, Wmax is the velocity at _Reff an d 
Q IS a power- law index. In Table 5 of ISteidel et"all (|20lO), 
the authors fit spectral line observations (of Lya, CIV, CII, 
Sill, SilV) using different parameter values ranging between: 
70 ^ -Rcff ^ 250 kpc, and 650 ^ Wmax ^ 820 km/s. Assum- 
ing Tmin — 1 kpc, they find that combinations of a within 
the range 1.15 ^ a ^ 1.95 together with another covariant 
parameter, can adequately reproduce observed ISM line pro- 
files. In the absence of a unique set of parameter values, we 
choose the following for one of our simulations (run RVWa, 
323}: r,nin = Ih-^ kpc, R^s = lOO/i"^ kpc, = 800 

km/s, and a = 1.15. Figure [T] (black curve in the left panel) 
shows the relevant outfiow velocity as a function of galacto- 
centric radius for this choice of parameters. 

Such a trend of Vw{r) correlating positively with r 
has also been found in other observations: for example 
IVeilleux et al.l (|l994l ) parameterized outflow velocity as 
Vw OC r" with 2 < n < 3 for the spiral galaxy NGC3079. 

We note that a scenario in which the wind velocity is not 
unique but has instead any value (fast and/or slow outflow) 
could possibly result in faster ejecta at large distances from 
the galaxy center and slower ones closer to the center, since 
the former would have had more time to travel out to large 
distances. The model presented here could mimic such a 
behaviour. 

We explore a further halo mass dependent parametriza- 
tion of our new wind model in run RVWb f i]2.5|l . Obser- 
vational studies of ultr a lumin ous infrared galaxies at z — 
0.042 - 0.16 bv iMartinI |200i) detect a positive correlation 
of galactic outflow speed with galaxy mass. Their Figure 
7 (bottom panel) shows that the outflow terminal veloci- 
ties are always 2 — 3 times larger than the galactic rotation 
speed. Inspired by these results, we assume that the velocity 
Vmax is related to the circular velocity Vdic of the halo at the 
virial radius i?2oo as: 



— 2'ycirc — 2 



Gi\fhalo 
R2OO 



(6) 



The halo mass Mhaio and radius are related such that ii200 
encloses a density 200 times the mean density of the Universe 
at the simulation redshift z: 
4n . 



Mh,io = ^i?200 (200pcritfiA/,o) (1 + 



(7) 
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Figure 1. Left panel: Black curve is the radial profile o f outflow veloc ity as formulated in Eq. (O, which is implemented in our 
hydrodynamic simulations motivated by observations of Stei del et al ] l|201Ch . The parameter values are for run RVWa (Table fljl : rmin — 

^max — 800 km/s, and a — 1.15. The blue horizontal straight line denotes the constant value of velocity, 
V = 400 km/s, in run CW. Right panel: Velocity at i?cfl as a function of halo mass, shown in E qs. l[6|l and |[8]l . for different redshifts of 
collapse. Such a formulation is implemented in run RVWb motivated by observational studies of iMartinI | |2005| ). 



Here pcrit = 3Ho / (SivG) is the present critical density. The 
assumption during code implementation is that the total 
group mass obtained by our group finder (i ]2.4p is the virial 
halo mass and that it is collapsing at the simulation red- 
shift. Using such a prescription the resulting dependence of 
maximum velocity on halo mass and redshift is: 

v.^..<xMlilil + z)'^\ (8) 

We also set the effective radius as: 

Rcs = i?200 oc M^il {l + zr\ (9) 

because Umax is the model wind velocity at Rch . Figure [T] 
(right panel) shows this parameterization of Umax as a func- 
tion of halo mass for different collapse redshifts. 

2.4 Implementation in the GADGET-3 code 

Our new wind model described in H2.3l reauires an extra pa- 
rameter: the distance of gas particles from their host galaxy 
center. We identify galaxies by running a Friends-of-Friends 
(FOF) group finder on-the-fly within our GADGET-3 simu- 
lations. Given a linking length, Lgr, the FOF algorithm finds 
a group by selecting all particles that lie within a distance 
Lgr from any other particle in the group. 

Firstly, an average inter-particle distance, Lp, is cal- 
culated; using the total mass, number of particles and the 
mean density. The linking length for the default on-the-fly 
FOF in GADGET-3 is Lgr/Lp = 0.16, which is commonly 
used to find dark matter (DM) halos. However, sometimes 
the FOF can link together DM substructures which actually 
belong to multiple galaxies with distinct stellar structures, 
into a single halo. Therefore, the stellar components of the 
FOF groups are considered to be more stable and thereby 
identified as galaxies. 

In order to select galaxies, we obtain stellar groups by 
modifying the FOF group finder as follows: we allow the 
FOF to first link over stars as the primary particle type, 
and then link over (gas + DM) particles as the secondary 
type. To facilitate selection of stellar groups, which are gen- 
erally smaller in size, we use a linking length three times 
smaller than the default one: Lgr/Lp — 0.16/3. The FOF 



is executed on-the-fly within a simulation at time intervals 
of a multiplicative factor 1.001 of the scale factor a, or, 
ctnext/ciprev = 1.001. All the groups abovc a minimum length 
of ^ 32 particles are finally selected, which provide us with 
a set of galaxies at a given simulation time. 

Our preliminary tests indicate that the center-of-mass 
position of the FOF groups are offset from the gas density 
peak location. Hence we decide not to consider the center- 
of-mass positions but to record the position of the member 
gas particle with maximum SPH density in each group: this 
maximum gas density location is considered as the group 
center. We tag in the code the group (representing a galaxy) 
and the group center each relevant gas particle belongs to. 
We define the galactocentric radius, r, as the distance of 
each gas particle from its group center. 

Assuming that outflows are generated from mass loss 
during SF, the implementation of launching the wind in- 
volves a probabilistic criteria similar to spawning star par- 
ticles, as described in SH03. Gas particles undergoing SF 
are stochastically selected and kicked into wind, by impart- 
ing an one-time boost. A probability is calculated in a 
timestep At for each multiphase gas particle: 

p=l-exp(-!|^). (10) 

Here x is the mass fraction of cold clouds and tgfr is the 
SF timescale, which depends on the parameters of the SF 
model. A random number is drawn in the interval [0, 1], and 
if it falls below p then the particle is given a wind kick. If 
V is particle velocity and (j> its gravitational potential, then 
its velocity after receiving wind kick is updated to: 

Uncw = «old + V^y- (11) 

The direction (unit vector y) is set as a random orientation 
along Void X V(^, which makes the wind particles to be pref- 
erentially ejected along the rotation axis of the galaxy or 
perpendicular to the galaxy disk. 

In order to enable the outflow to escape from dense 
SF regions without affecting the SF, a new wind particle 
(gas particle just receiving a wind kick following Eq. Ilip is 
decoupled from hydrodynamic interactions, for a maximum 
duration tdec = 0.025tH(2), where tu{z) is the Hubble time 
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at a relevant simulation redshift. Then the wind particle does 
not enter hydrodynamic force computation in the code, but 
is included in gravity and SPH density calculation. If the 
particle's density has fallen below pdcc = 0.25/9sf, then the 
decoupling is stopped before tdcc, and full hydrodynamics 
are enabled again. 

Summarizing, the main free parameters in the wind 
model are and r), along with tdoc and pdcc for the de- 
coupling prescription. The wind velocity is either constant 
(Eq. O i|2.2|l , or dependent on radius (Eq. [5]) and halo mass 
(Eq.[8j. 

In this Section we have shown the actual implemen- 
tation of different forms of Vw, including radial and mass 
dependence of velocity which are observationally motivated. 
We later study the implications such models have on the 
CGM and IGM properties in the simulated volume. 



2.5 Simulations 

The series of simulations we perform are listed in Table [T] A 
cosmological volume, with periodic boundary conditions, is 
evolved starting with an equal number of DM and gas parti- 
cles from 2; = 99, where the initial co nditions have been gen- 
erated using the CAMlfl software l|Lewis et al.l bOOOD . up 
to z = 2. A concordance flat ACDM model is used with the 
following parameters: Q.Mfi — 0.2711, Qa,o ~ 0.7289, Qb,o ~ 
0.0463, ns = 0.96, as = 0.809, i/o = 70.3 km s*^ Mpc"\ 
in agreement with recent observations of the cosmic mi- 
crowave background radiation, weak gravitational lensing, 
Lyman- Q f orest and mass function evolution of galaxy clus- 
ters (e.g.. iLesgourgues et al.l l2007l : IVikhlinin et al. I l2009l : 
iKomatsu et al.ll201ll '). 

The smaller box-size Lbox ~ 5h~^ Mpc comoving series 
(or, box SB) has A''part = 2 x 128^ DM and gas particles 
in the initial condition. This is the higher resolution series 
with gas particle mass mgas = 7.66 X W^h'^MQ, and DM 
particle mass moM = 4.49 X lO^h-^MQ. The Plummer- 
equivalent softening length for gravitational forces is set 
to -Lsott ~ 0.98/;."^ kpc comoving. In the larger box-size 
ibox = 25/i~^ Mpc comoving series (or, box LB) the num- 



bers are: iVpart = 320'', mgas = 6.13 x lO'^/i'^M©, moA/ = 
3.59 x lO''/i"^M0,Lsoft = 1.95/i"^ kpc. The minimum gas 
smoothing length is set to a fraction 0.001 of L^oit in all the 
simulations. The second set of simulations, at slightly worse 
resolution compared to the first set, has been chosen in or- 
der to increase the statistics, since more number of halos 
and more massive halos are expected to form in the larger 
box. 

In each series four runs are performed incorporating the 
same non-wind subgrid physics described in >i2.2l and inves- 
tigating different galactic wind models: 

(1) NW: no wind; 

(2) CW: energy-driven wind with constant velocity ( H2.2I 
Eq. O t;„ = 400 km/s; 

(3) RVWa: radially varying wind with fixed parameters 
( ^ r„»„ = Ih-^ kpc, R,ff = imh~^ kpc, = 800 
km/s, a — 1.15; and 

(4) RVWb: radially varying wind with parameters depen- 
dent on halo mass. 



http: / /camb.info/ 



The smaller box-size runs have an extra t appended to the 
names. Few results are presented for both the box sizes, 
namely SB and LB runs. While the remaining results are 
shown only for one box-size, with the reason given in the 
relevant parts of 33] next. 

The new velocity prescriptions are plotted in Figure [T] 
Left panel shows the outfiow velocity radial profile in run 
RVWa as the black curve, and for comparison, the blue hor- 
izontal line denotes the constant value of velocity, v — 400 
km/s, used in run CW. The right panel depicts the velocity 
at i?cff in run RVWb as a function of halo mass for different 
redshifts of collapse. 

We have analyzed the carbon content of the gas in the 
box, or metallicity of the element carbon, Zc- It is computed 
as the ratio of carbon mass to the total particle mass for 
each gas particle. Abundance ratios are expressed in terms 
of the Solar metallicity, which is Zc,q = 0.002177 (mass 
fraction of carbon in Sun) derived from the compilation by 
lAsplund. Grevesse fc Sauva 1 (|2005l ). 



3 RESULTS AND DISCUSSION 
3.1 Outflow Velocity 

The different implementations of the outflow velocity in 
the SB runs are shown in the panels of Fig. (2] NWt (top- 
left), CWt (top-right), RVWat (bottom- left), and RVWbt 
(bottom-right). Here the smaller boxsize is chosen for visu- 
alization purposes, to unambiguously distinguish the wind 
velocity from the remaining gas in the most-massive galaxy 
within the box volume which has the largest SFR. The ve- 
locity magnitude of the gas particles is plotted as a function 
of distance from center of galaxy at 2 = 2.44. As described 
in i|2.4l galaxies are identified on-the-fiy in the simulations 
using a FOF algorithm linking over star particles as the pri- 
mary type using a linking length 3 times smaller than the 
default one, and the galaxy center corresponds to the loca- 
tion of the maximum gas density. 

All the gas particles lying within a radius of WOh~^ kpc 
from group center are shown as black points in Fig. [21 The 
particles which have recently received a wind kick velocity 
can be tracked, for the duration of their hydrodynamic de- 
coupling, through their positive decoupling time ( tj2.4p . We 
call these the wind-phase particles and we mark them with 
the red plus symbols. The median wind-phase particle's ve- 
locity in radial bins is plotted as the green asterisks joined 
by the green dashed curve. The input subgrid wind velocity 
(vw) implemented in the code is instead represented by the 
blue solid curve. In order to mark the differences clearly, in 
each run results are stacked over the 20 most-massive galax- 
ies having total halo (DM -I- gas + stars) mass between 
A/halo = (0.9 — 15) X 10^" M0, the exact range for each case 
is written in the respective panel. 

Run NWt gives the radial structure of gas velocity 
without any wind, reaching a maximum of u ~ 400 — 500 
km/s around galaxies. The other 3 runs have winds im- 
plemented and present a distribution of red points. Wind 
particles are kicked over a radial range from the center up 
to r ~ (30 — 50)h~^ kpc, marking the size of star-forming 
galaxy disk. These particles gain an additional velocity of 
Viu during their wind kick ( ^2.4p . Runs CWt and RVWat 
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Figure 2. Velocity magnitude of the gas as a function of distance from galaxy center at 2 = 2.44 in the SB runs with different wind 
models: no wind (top- left), energy-driven constant- velocity wind (top-right), radially varying wind with fixed parameters (bottom- left), 
and radially varying wind with parameters dependent on halo mass (bottom-right). All the gas particles within 100/i~^ kpc from the 
group center are shown as black points. In the runs with wind, the red plus symbols are the wind-phase particles (described in ^3.ip . 
whose median velocity in radial bins is denoted by the green asterisks joined by the green dashed curve. The blue solid curve represents 
the subgrid wind velocity implemented in the code in each case; in the last panel (run RVWbt) the blue solid line refers to an intermediate 
halo mass Mj^aio = 8 X 10^^ Mq. Each panel shows the stacked result of the 20 most-massive galaxies having total mass within the range 
specified in the panels. 



demonstrate this as some red points clustered around or 
somewhat above the blue curve. Afterward they move away 
from galaxy center and get decelerated, causing a fraction 
of the red points to fall well below Vw Therefore, when con- 
sidering the wind particles at given r, the median veloc- 
ity (green dashed curve) is somewhat lower than (blue). 
Run CWt has particles kicked by = 400 km/s at all 
r, while in our new implementation RVWat the upper en- 
velope of the red points follow the Vw{r) radial profile of 
Eq. [S] In run RVWbt (bottom-right panel), the blue solid 
curve shows Vw for a representative intermediate halo mass 
Mhaio = 8 X IO^^Mq. This turns out to be a weak wind, 
because most of the gas is moving at velocity above the blue 
curve at r 10/i~^ kpc. Hence the wind particles do not 
move outward appreciably, and end up following a distribu- 
tion very close to the remaining gas (black points). 

Fig. [2] also demonstrates that even though wind par- 
ticles are kicked only one time according to the original 
prescription of SH03, since this procedure is done at every 
timestep (stochastically selecting from all the star-forming 
particles), the desired is reproduced by multiple gas par- 
ticles at any given time. The wind particles slow down, but 
in the next timesteps more particles are kicked out with Vw, 



as long as there is SF, thus roughly mimicking a continu- 
ous outflow of gas moving with given subgrid input speed 
around a galaxy. 

Velocity magnitude histograms of gas particles at 2: = 
2.91 are plotted in Fig. |31 showing mass fraction in velocity 
bins; wind (gas) particles are represented by solid (dashed) 
curves. In each case the histogram is normalized to the to- 
tal gas mass in the respective simulation. The SB runs are 
shown in the left panel where the maximum particle veloc- 
ity is ~ 700 km/s, and LB runs in the right panel where 
particles reach ~ 1100 km/s because of higher-mass halos 
forming in a larger box. 

Runs CWt and CW have peaks in their wind particle 
velocity histogram at w ~ 400 km/s, since in this case Vw is 
kept constant at this value. The wind velocity histogram is 
much flatter in runs RVWat and RVWa with no well-defined 
peak, because has a range of values depending on galacto- 
centric distance. Runs RVWbt and RVWb also have a peak 
wind velocity, which is at a lower value v ~ 100 — 200 km/s, 
corresponding to the outfiow of lower-mass halos which are 
more numerous in hierarchical structure formation. Consid- 
ering all the gas particles, the no- wind (NW) run has slightly 
higher mass fraction at ti ^ 300 km/s (not distinguishable in 
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Figure 3. Velocity magnitude histogram of wind particles (solid), and all gas particles (dashed) at 2: = 2.91, showing the mass fraction 
of gas per velocity bins. The SB runs are in the left panel and LB runs are in the right panel, in each case the histogram is normalized 
to the total gas mass in the simulation at the redshift considered. 



the plot), while runs with wind (CW, RVWa, RVWb) have 
higher gas mass with v ^ 300 km/s. 



3.2 Global Star Formation Rate Density 

The global star formation rate density (SFRD) as a function 
of redshift is plotted in Fig. 14] for the SB runs in the left 
panel, and for the LB runs in the right one; in each panel 
the four wind models are labeled by the different colors and 
plotting symbols. The SFRD is computed by summing over 
all the SF occurring in the whole simulation box and dividing 
it by the time-step interval and box volume to obtain the 
rate density in Mq yr~^ Mpc~^. Galactic wind feedback 
clearly has significant impact on SFRD, reducing SF several 
times depending on the outflow model. 

Between z ~ 8 — 10, there is already a trend of sup- 
pression of SFRD with winds as the CW run has 1.2 times 
lower SFRD than NW, while at z < 8 the SFRDs of the 
runs diverge more. The differences are small between runs 
RVWb and NW; in RVWb the SFRD values are 1.5 - 2 
times smaller than NW between z « 3 — 8, and agree with 
NW at 2; « 2 — 3, making RVWb the wind model least ef- 
fective in suppressing SF. There is however a substantial 
reduction of SFRD in runs CW and RVWa. The SB box 
have a stronger wind-driven suppression; SFRD is 12 times 
smaller in RVWat (and 5 times smaller in CWt) than NWt 
at 2 < 3. While in the LB box, the reduction is 4 times in 
RVWa and 2 times in CW at 2 < 4. Analyzing the effective 
wind models, RVWa has 1.2 times higher SFRD than CW 
from high z up to z = 5, but later RVWa suppresses it more 
and produces 2 — 4 times lower SFRD aX z = 2. 

Effects of box-size and resolution are visible in our re- 
sults. The z = 12 SFRD is ~ 2 higher in the SB box than 
LB, because the SB series has higher resolution, and can 
track denser gas, forming more stars. However at 2 < 3 — 4, 
the RVWat run of SB series produces 2 — 3 times lower 
SFRD than RVWa of LB, one reason for which is the rele- 
vant box-size; the halos forming in box SB are less-massive 
and fewer in number than in LB. Studying the shape of the 
SFRD evolution, the SB runs tend to show a peak in SFRD 



at a certain redshift: z = 2.5 in NWt and RVWbt, 2 ^ 2 in 
CWt, and z = 4 in RVWat. Whereas all the LB runs have 
a plateau of SFRD between z = 2 — 3.5. 

Observational data are overplotted in Fig. |3] with 
the cyan symbols and error bars. Each data set is 
shown with a different plotting s ymbol as listed next. 
These data are taken mainly from ICucciati et al.l (|2012| ) 
- filled ci rcles, and t h e co mpilations th erein origi- 
nally from ISteidel et all (Il999l) - asterisks, lOuchi et al.l 



(|2004h - plus si gns. iPerez-Gonzalez et al 
verted triangle s, Schiminovich et al.l (|2005 )^ 



20051 ) - m- 

- diamonds. 



ISchiminovich et al.l 

iBouwens etal] (|2009l) - open squares, iReddv fc Steidell 
||2009D - crosses, |R .odighier o et al.l ([2010') - open circles, 
Ivan der Burg. Hildebrandt fc ErberJ ( 2010 ) - upright tnan- 
ote5, iBouwens et al.l l|2012l ) - filled squares. 

Comparing with observations, we see that no single sim- 
ulation model can fit the data from a single observation. 
Taken collectively, there is a better match of simulations 
with observations at low-z. At earlier cosmic epochs between 
z ~ 4.5 — 10, SFRD in the simulations is systematically 
higher, reaching 2 — 10 times the observed values. At later 
times between z ^ 2—4.5, most of the observations lie within 
the ranges of SFRD produced by the different wind models. 
Run CW of the LB series is the model providing the best-fit 
to the observations at lower-z. 



3.3 Galaxy Mass Function and Mass Fraction 

The gas and stellar mass function {dN, in units of (log M)~^ 
Mpc~^) of galaxies in the LB runs at 2 = 2.23 is plotted as 
the solid curves in Fig. [5] left panels. The occurrence of a 
peak of dN at an intermediate mass, and the drop of dN at 
subsequent lower masses is due to resolution effects. The top- 
left panel shows the mass function of the gas component in 
galaxies. All the 4 runs have the same slope between Mgas 
(10^ — 1O^^)M0, but shifted sideways as described next. NW 
and RVWb produce almost similar trends: dN « 0.3 at (2 — 
3) X 10^ Mq, with the RVWb mass function shifted right- 
ward by ~ 1.2 mass units (Mq). Results in CW and RVWa 
are quite similar, shifted left-ward with respect to NW by 
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Figure 4. Star formation rate density in whole simulation volume as a function of redshift, for the SB runs in the left and LB in the right 
panel, wit h the respective wind models labeled by the color and plotting symbol. The cyan symbols and error b ars denote observational 
data from ICucciati et alJ (l2012h - filled c i rcles, and the compilation s therein originally f r omlSteidel et aD lll999^ - asterisks , [Ouc hi et alj 
plus signs. ^Perez-Gon zalez et ahl | |2005|) - inverted tr i angles , ISchiminovich e t al ? ('2005') - diamonds, 'Bouwens e t al"] ll2009t) - open 
squares, |Rcddv & Stcidcl ( 2003) - crosses, iRodighiero et al - open circles, .van dcr Burg. Hildcbrandt fc Erbcnl l|2010r i - upright 

triangles, I Bouwens et al.i ||2012| ) - filled squares. Detailed comparison is in i|3.2l 



(2 — 3)M0. The wind in these runs expel gas efficiently out of 
galaxies with Mgas 55 10^ Mq, producing a smaller number of 
objects with such high gas masses. Only at J5 10^'^ Mq, CW 
has slightly more galaxies at the same mass than RVWa. 

The bottom-left panel of Fig. [5] shows the mass func- 
tion of the galactic stellar component. The sharp cutoff at 
Mstciiar « 5 X 1O^M0 corresponds to the minimum mass 
of stellar groups found by the FOF group finder algorithm 
in our simulations (32 star particles, i]2.4l) . Runs CW and 
RVWa having dA^ « 0.2-0.3 at Msteiiar = (1.5-2) x 10* M©, 
produce a steeper stellar mass function than NW and RVWb 
where diV f» 0.1 at lO^Afg. There is a small excess of galax- 
ies in RVWa compared to CW between Mstdiar = (2 x 10* — 
1O^)M0, and a slight excess in NW than RVWb in the range 
(10^ - 1O")M0. The stellar mass functions can be com- 
pared to semi-analytic modeling of lBower. Benson fc CrainI 
l|2012h . At Msteiiar = 10* Mq their Fig. 1 result at low-z 
lies between our CW / RVWa and NW / RVWb cases, at 
10^ Mq it is closer to our CW / RVWa, while at 10^° Mq is 
closer to our NW / RVWb results. Overall o ur mas s func - 
tions are steeper than that of lBower. Benson fc CrainI (|2012l ) 
and the obse rvational data therein, also other low-z obser- 
vations (e.g.. iPapastergis et al.|[2012l '). 

Fig.[5]also shows results of the wind model RVWat from 
the SB box, plotted as the red dashed curve in each panel. 
Comparing RVWa and RVWat, the smaller higher-resolution 
box has lower-mass halos forming, extending the gas mass 
function to A/gas ~ 10* Af©, and the stellar mass function 

down to A/stellar ~ 4 X IQ^Mq. 

The mass fractions of gas and stars in galaxies with 
respect to total mass of halos are plotted in the right panels 
of Fig. (5] gas mass fraction in the top-right, and stellar mass 
fraction in the bottom-right. Solid curves denote the median 
value within a mass bin for each run labeled by the color and 
plotting symbol. The grey shaded area encloses the 70th 
percentiles above and below the median in run RVWa (red 



curve), and show the typical scatter at given halo mass. 
With no-wind (NW) the gas fraction in galaxies decrease 
from A^gas/Afhaio ~ 0.13 at A4aio = IO^^A^q to 0.04 at 
4 X IQ^'^Mq. The wind (CW, RVWa, RVWb) flattens the 
Afgas/Aihaio trend making it oscillatory around (0.08 — 0.1) 
over the same Afhaio range. The gas fractions in galaxies are 
lower than the cosmic baryon to DM fraction used in our 
simulations {Q,b.o/^m,o ~ 0.17). The gas mass fractions we 
obtain are roughly consistent with those found in massive 
cluster halos (Afsoo > W^'^Mq) at 2 = 0, in simulations 
including baryonic ph ysics of cooling, SF and ACN feedback 
(jPlanelles et al.ll2012l '). 

The stellar fraction is largest in NW, because the SF 
rate is higher in the absence of wind ( H3.2|) producing more 
stars for the same halo mass. The wind in the other runs 
expel some star-forming gas out of galaxies, reducing the SF 
and consequently the stellar mass. A/gtciiar /-^^haio incrGcLses 
from « (0.01 - 0.025) at Afhaio = W^°Mq to « (0.08 - 0.16) 
at 6 X 10^^ Mq, at a varying rate for the different runs, and 
remain almost flat at higher halo masses up to 4 x 10^^ Mq. 
NW and RVWb have the steepest increase, followed by CW, 
and RVWa is the flattest. Observations also find the stellar 
mass fractions of \ow-z galaxies increasi ng with halo mass 
betw een Afhaio ~ (10^° - 10^^)AfQ (e.g., [Papastergis et al.1 
with a comparable slope, but the absolute values in 
observations are ~ 10 times smaller than our simulations. 

We thus conclude that the outflow models CW and 
RVWa have significant impact on the gas and stellar mass 
functions, as well as on the mass fractions of gas and stars 
in halos. 



3.4 Temperature-Density Phase Diagram 

The temperature versus density phase diagram of all the gas 
at 2: = 1.98 is plotted in Fig. (6] for the LB runs in 4 panels, 
color coded by the gas mass fraction. The larger boxsize 
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Figure 5. Galaxy gas mass function (top-left), stellar mass function (bottom- left), along with gas mass fraction (top-right) and stellar 
mass fraction (bottom-right) w.r.t. total mass of halos, at ^ = 2.23, of the LB runs shown as the solid curves. In the right panels the solid 
curves denote the median value within a mass bin for each run, and the grey shaded area enclose the 70 percentiles above and below the 
median in run RVWa (red curve) showing the typical scatter. One of the wind model (RVWat) from the SB box is plotted as the red 
dashed curve in each panel. Details are given in the text in ij3.3l 



is chosen for the two [T — p] diagrams we present in this 
section, because it gives a statistics over a larger range of 
halos. The a;— axis denotes the overdensity, 5 = Pgas/{pB), 
which is the comoving density of gas pgas as a fraction of the 
mean baryon density, {ps) = Qb,oSHq /{SttG). The vertical 
dashed hne is the SF threshold density (psF, i)2.1|) above 
which a fixed effective equation of state is imposed on the gas 
in the multiphase SF model (SH03, Wicrsma ct al. 2009). 
The hot-phase temperature is shown for the multiphase gas 
particles, which forms the constant-slope straight line above 
PSF in the top-right portion of each run. In the NW run (top- 
left panel), all the gas denser than psF follows the equation 
of state, and there is no cold gas at those high densities. 

In the runs with wind, a fraction of the gas denser than 
psF is cooler (< 10* K) than the NW case, forming a cold, 
dense tail in the bottom-right portion of the phase diagram. 
This is a consequence of the decoupling formalism of the 
galactic outflow implementation ( i]2.4p . The wind particles 
are allowed to move away from SF regions and no longer 
follow the SF effective equation of state. They are however 
still very dense, where the cooling rate is high and cooling 
time short. Therefore they cool very fast, soon after being 



kicked into wind, and form the cold, dense extension above 
pSF in the phase diagram. Subsequently, if the kick veloc- 
ity is higher than the neighbors, the wind particles exit the 
dense SF region soon, moving to lower-density (< psf), hot- 
ter (> 10* - lO'^ K) regions. Thus in our galactic outflow 
models, the wind is launched from dense SF phase, and soon 
after that the wind goes through a cold phase. 

The remaining of the phase diagram looks qualitatively 
similar in the 4 runs. The majority of the gas (red in Fig. [6| 
lies in underdense to moderately overdense (5 ~ 10~^ — 10^) 
and cool-warm (T ~ 10"^ — 10^ K) phase, which is the cos- 
mic baryons in balance between adiabatic cooling and pho- 
toionization background heating. A small fraction of gas at 
a higher range of overdensities {S ~ 10~^ — 10*) is hot 
(r ~ lO'' — 10* K), being heated by shocks during gravi- 
tational collapse. At (5 ~ 10^ — 10*, there is another cold 
(r < 10''' K) phase which is the dense gas cooling into DM 
halos to form galaxies. 

The gas mass fraction in the phases varies between the 4 
runs, forming different patterns in the phase diagram. Runs 
NW and RVWb have a higher fraction of gas in the multi- 
phase SF branch of effective equation of state. Among the 
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Figure 6. Temperature versus overdensity (gas comoving density as a fraction of the mean baryon density) phase diagram of all the gas 
at 2 = 1.98, for the LB runs with name labeled in each panel. The color code represents the mass fraction (in log-scale) of gas in density 
- temperature bins, with red showing the highest fractions going to black for the lowest. The vertical dashed line is the star formation 
threshold density (physical psp fsi 0.1 cm~^, i|2.1ll above which a fixed equation of state is imposed on the gas in the multiphase SF 
model. The hot-phase temperature is shown for the multiphase gas particles, which results in the constant-slope straight line above psF 
in the top-right portion of each panel (" ^13.411 . 



wind runs, RVWb has the highest mass fraction in the cold 
wind phase, followed by CW and finally RVWa. This is be- 
cause of the different wind kick velocity in each case; RVWb 
has the lowest (seen in Fig. (2] i|3.1[l where gas is not able 
to escape the dense phase, while in RVWa the wind particles 
are able to escape the dense phase quickly. Runs CW and 
RVWa have a higher fraction of underdense (5 < 1), warm- 
hot (T ~ 10* — 10® K) gas, denoted by the red and yellow 
contours in Fig. [6l composed of gas previously ejected from 
dense SF regions of galaxies as strong winds. 

The temperature - density phase diagram of only the 
CGM gas aX z = 1.98 is plotted in Fig. [7] for the LB runs 
in 4 panels. Gas particles lying inside the virial radius ii200 
(computed analytically using Eq. [7|) of all the FOF halos are 
selected, excluding the star-forming ones (ie. particles with 
SFR = are only counted). The selected CGM gas comprise 
of a number fraction (0.07 - 0.1) of all the gas particles in 
a run, and is shown in Fig. [T] color coded by the gas mass 
in [T — S\ bins as a fraction of the total mass plotted. The 
densities defining the boundaries of the CGM are denoted 
by the vertical dashed lines. The SF threshold density is 
the upper limit, above which gas forms stars and becomes 
part of galaxy. From Fig. [71 we infer a lower limiting density 
fcoM = 5 of the CGM. 

The NW run (top-left panel) has no CGM gas denser 
than PSF, because of imposing the SFR = criterion for 
selecting CGM particles, as here all the gas denser than psp 
is star-forming. The runs with wind have a fraction of gas 
denser than psF which is not star-forming, and they are the 
wind-phase particles: mostly cold undergoing hydrodynamic 



decoupling, and a small fraction which recently received a 
wind kick and still in the hot multiphase branch. Noting the 
effect of feedback, we see that the CGM is separated into a 
cold T < 10^ K phase and hotter gas in a more well-defined 
way in models CW and RVWa impacted by their effective 
outfiows, than in NW and RVWb. 

Overall, analyzing the [T — S\ plane, we can conclude 
that the thermal properties of the IGM remain roughly the 
same between the different outfiow models; with the winds 
producing a higher fraction of underdense, warm-hot gas. 
Furthermore the CGM is well separated into cold and warm- 
hot phases by the winds. 

3.5 Single Galaxy Gas Kinematics 

The most-massive galaxies at z = 2.12 in the SB simula- 
tions are plotted in Fig. [8] The smaller boxsize is used here 
again for ease of visualization (like in ^3.ip . to clearly dis- 
tinguish the outfiowing gas in the plot of a massive galaxy 
within the box volume. Each of the 4 rows shows a unique 
gas property, the different wind model runs are labelled in 
the top-left corner of the 1st, 3rd and 4th rows. The max- 
imum gas density location is at the center, and the pan- 
els show the projection of gas kinematics in the [x — y\ 
plane of a surrounding {2QQh~^ kpc)"^ volume. The total 
halo mass (DM -f gas + stars) of the galaxies are in the 
range Mhaio = (1.4— 1.9) x 10^^ Mq. One can see the forma- 
tion of a galaxy disk in runs NWt, RVWat and RVWbt of 
varying sizes, whereas in CWt the galaxy looks irregularly 
shaped in the projected plane. 
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Figure 7. Temperature versus overdensity phase diagram of tlie CGM gas at z = 1.98, for ttie LB runs, in a format similar to Fig.|6] 
The gas representing the CGM is selected by tracking particles lying inside -R200 of all the halos and excluding the star-forming ones, 
which comprise of a number fraction (0.07 - 0.1) of all the gas particles in a run (all gas was shown in Fig. [SJ. The vertical dashed 
lines represent the bounding densities of the CGM: the star formation threshold density is the upper limit, and a lower limiting density 
5cGM = 5 inferred from this plot. 



In the top two rows the outflowing (radial velocity with 
respect to galaxy center, Vr > 0) particles are denoted as 
red, and the inflowing (vr < 0) as black. The first row shows 
the positions of all the gas particles within the projected 
volume, and the second row depicts the velocity vectors of 
10% outflowing gas particles and 10% inflowing gas. The 
number fraction of outflowing (/out) and inflowing (/in) gas 
are written in the top-left corner of the second row panels. 
Within this {200h~^ kpc)^ projected volume, in runs NWt, 
RVWat and RVWbt, about half of the gas is inflowing and 
the other half outflowing. While in run CWt more gas (64%) 
is undergoing infall than outflow (3 6%). These f ractio ns in 
run CWt are similar to the result bv lShen et aP (|2012h who 
found that in the Eris2 zoom-in simulation, the CGM of a 
Milky Way-type galaxy at z ~ 3 has about one third of all 
the gas within the virial radius as outflowing. In our Fig. |8l 
the inflow is overplotted on the outflow in the flrst row de- 
picting a prominence of black points over red at the center, 
and vice versa in the second row showing a prominence of 
red points over black; this demonstrates the mixture of gas 
dynamics (i.e. contains both inflow and outflow) in the cen- 
tral 20 — 30 kpc regions. 

The outflow is not well-structured in the no-wind 
(NWt) and the low-velocity wind model (RVWbt) cases, 
some gas escaping along the direction of least resistance 
through a low-density void in the bottom-right half. In run 
CWt most of the outflowing gas lies inside r < 50 kpc, be- 
cause here the wind kick velocity Vw = 400 km/s is not 
enough to drive the outflow to larger distances; the strong 
gravitational potential of the halo causes the gas to fall back 
inward also resulting in a higher inflow fraction in this run. 



The radially- varying subgrid wind velocity v^{r) of Eq. [5] 
in run RVWat produces a well-developed gas outflow propa- 
gating perpendicular to the galaxy disk, escaping to r > 100 
kpc from the center, seen as the red arrows along top-right 
and bottom-left in the 3rd panel of the 2nd row. 

The third row of Fig. [8] shows the gas overdensity color 
coded on a log-scale from red as the highest and black as 
the lowest values. All the runs have a central overdense re- 
gion, where SF occurs forming a galaxy disk in some cases. 
Wind particles are kicked from these SF regions, carrying a 
fraction of the central gas to larger distances. In run CWt 
the constant-velocity wind forms an extended halo of gas at 
r ^ 50 kpc, seen as yellow-green in the 2nd panel, which is 
most likely bound to the galaxy potential and not able to es- 
cape. The strong radially-varying wind in RVWat produces 
a more extended (r ^ 100 kpc) but lower-density diffuse 
outflow, which is likely escaping the galaxy in the 3rd panel. 
The low-velocity wind run RVWbt has most of the gas con- 
centrated in the central 10 — 20 kpc. 

Carbon metallicity {Zc, deflned in i )2.5p is plotted in the 
fourth row of Fig. [S] Runs NWt and RVWbt have a larger 
central concentration of metals (red in the flgure) originating 
from SF, because there is either no-wind or the wind is not 
effective to spread the metals around. The central metal 
content is smaller in runs CWt and RVWat, since winds 
carry the metals out from the SF regions and enrich the 
CGM and IGM. Consequently, in the 2nd and 3rd panels, 
the distribution of metals outside galaxy disk is aligned with 
the outflowing velocity (red arrows in the 2nd row) , and two 
extended gas outflow regions (overdensity in the 3rd row) are 
present. 
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Figure 8. Projection of gas kinematics in tlie {x — y] plane of a (200/i~^ kpc)^ volume centered around the maximum gas density 
location of the most-massive galaxy at redshift z = 2.12, in the SB simulations with different wind models (from the left): no wind, 
energy-driven constant-velocity, radially varying with fixed parameters, and that with halo mass dependent parameters. In the top two 
rows the outflowing {vr > 0) particles are denoted as red, and the inflowing {vr < 0) as black. First row from the top shows the positions 
of all the gas particles within the projected volume, and the second row depicts the velocity vectors of 10% of each component. The total 
halo mass of each galaxy, and the number fraction of outflowing (/out) and inflowing (/in) gas are written within the panels. The inflow 
is overplotted on the outflow in the first row, and vice versa in the second. Third row shows gas overdensity, and fourth row is carbon 
metallicity, color coded from red as the highest and black as the lowest values. 
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The large scale structure of all the runs are similar as 
expected from the same initial condition, though there are 
small-scale differences as described above arising from differ- 
ent wind models. There is gas inflow from bottom / bottom- 
right regions of the panels into the center and from top / 
top-right into the center, seen as black points in the 1st row, 
black inflowing arrows in the 2nd, overdense fllaments in the 
3rd, which have lower-metallicity as seen in the bottom row. 
These likely denote pristine gas infall along cosmological flla- 
ments that "feed" the galaxy at the center. Another common 
feature is the overdense gas clump 70 — 80 kpc to the top of 
the center, most likely an infalling sub-structure which will 
eventually merge with the central galaxy. 



3.6 Radial Profiles Around Galaxy Centers at 

2 ~ 2 

The radial proflles of gas properties around centers of galax- 
ies (found by FOF group flnder) a,t z = 1.98 for the LB runs 
are plotted in the next set of figures (Figs.[ni[Tni[IIl[IlJ. The 
larger boxsize is chosen for all the radial profiles we present 
in this section, since more number of halos and more massive 
halos form in the larger box; it hence gives a wider statistics. 

Each of these figures show a property as a function of 
comoving radius (or distance from the maximum density 
position considered as the galaxy center), by counting all 
the gas particles lying inside 300 kpc from the center, 
in the following format. The four panels denote total (DM 
-|- gas -I- star) halo mass (Mhaio/A^©) ranges: 10^ — 10^° 
(top-left) with number of halos in the four plotted runs 
within the range iVhaio = 739 - 997; 10^° - lO" (top- 
right) having iVhaio = 7453 - 7906; lO" - lO^^ (bottom- 
left) where iVhaio = 270 - 393; 10^^ - lO" (bottom-right) 
with A^haio = 10 — 12. All the halos within each mass range 
are stacked, and the plotted solid curve denotes the median 
value in radial bins for each run labeled by the color and 
plotting symbol. The grey shaded area encloses the 70th 
percentiles above and below the median in run RVWa (red 
curve), showing the typical radial scatter at a given radius, 
since galaxies in general do not have spherically-symmetric 
properties. The vertical dashed line is the halo virial radius 
-R200 in comoving coordinates (analytical expression from 
Eq. [7]) for the following masses in the panels from top- left: 
Mhaio/Me = 3 X 10^ 3 x 10^°, 3 x lO" and 3 x 10l^ where 
the exact values are R200 = 32.2, 69.4, 149.6, and 322.2 
kpc respectively. 

3. 6. 1 Density 

The gas overdensity (ratio of density to the mean baryon 
density, ^3.4p radial profiles are plotted in Fig. [9] Within 
the approximate virial radius r < R200 all the proflles are 
composed of two negative-sloped functions separated by a 
threshold radius dependent on halo mass and wind model. 
At r > -R200 the density proflles in the two lower halo mass 
ranges (top two panels) rise again and fall, forming a local 
peak. This occurs because of the presence of other smaller 
halos and substructures, giving rise to a local density peak. 
This trend is visible only around lower mass halos since for 
those a sufficiently large volume within r ~ (4 — 9)-R2oo is 
plotted to reach surrounding structures. 



The inner parts r < 10h~ kpc of the 2 lower halo mass 
ranges (10^ — 10^*^ and 10^*^ — 10^^, top two panels) present 
significant differences: CW and RVWa produce a lower den- 
sity than the NW and RVWb cases, because the strong wind 
is able to expel gas from the star-forming regions, conse- 
quently reducing the central gas density by 10 — 30 times. 
The trend is almost reversed in the outer parts (r > -R200): 
CW, RVWa and RVWb has a somewhat higher density than 
NW, because of accumulation of gas expelled by wind. 

The differences are smaller in the two higher halo mass 
ranges (lO" - 10^^ and 10^^ - 10^^, bottom two panels), be- 
cause the wind, not being dependent on halo mass in runs 
CW and RVWa, is less effective in ejecting gas out of mas- 
sive galaxies. CW produces slightly higher density (by 2 — 3 
times) than RVWa in the inner r < 10h~^ kpc, while outside 
this distance they are quite similar. 



3. 6. 2 Temperature 

The temperature radial profiles are presented in Fig. IIOI 
where the hot-phase temperature has been used for those 
gas particles which are multiphase (star- forming). In the in- 
ner parts r ^ 10/i~^ kpc of the galaxies, the T-profiles follow 
the negative-sloped density-profiles (Fig. [9]), for run NW in 
all the panels, run CW in top-right and bottom-left, and 
run RVWb in top-right. This represents the dense gas near 
galaxy center undergoing SF, and having a high tempera- 
ture (~ 10^ — 10^ K) as a result of following the SF effective 
equation of state ( H3.4p . In the wind runs there is also the 
presence of cooler (< 10* K) gas, which has recently received 
a wind kick and undergoes hydrodynamic decoupling ( i]2.4p . 
This gas was seen as the cold, dense tail in the bottom-right 
portion of the [T — p] phase diagram (Fig. [6l i|3.4p . Conse- 
quently, runs CW, RVWa and RVWb present a bimodal T 
distribution at r ^ 10h~^ kpc in most of the panels, com- 
posed of hot multiphase SF and cold wind. 

There is a change in T slope in the outer parts, at r ^ 
6h~^ kpc in the top 2 panels and at r ^ (20 — 30)/i~^ kpc in 
the bottom 2 panels, when the gas T increases with radius, 
likely because of shock heating at galaxy outskirts. Between 
(10 - 300)/i"^ kpc, the NW model results in the highest-T 
(~ lO''' - lO'' K), followed by RVWb, then RVWa and CW. 

3.6.3 Carbon Metallicity 

Radial proflles of carbon metallicity are plotted in Fig. Illl 
showing the ratio of carbon mass fraction in the gas to that 
of the Sun. The solid curve medians and percentile values are 
computed considering all (both enriched and non-enriched) 
gas particles in radial bins. The dashed curves represent 
median- for the enriched particles only, i.e. those having 
Zc > 0, without counting particles with Zc = 0. 

We also checked the gas fraction (by particle number) 
which has Zc > 0. All the gas is enriched inside a limiting 
radius, which is ~ 7, 10,40, and 120 kpc in the panels 
from top- left, deflning the size of the central SF region or 
the sites of metal generation. The enriched fraction decreases 
outside the limiting radius, with a trend depending on halo 
mass and wind model: it reduces to (0.5 — 0.8) at -R200 and 
(0.45 - 0.7) at 300 h'^ kpc. Runs CW and RVWa enrich a 
higher (up to 1.4 times) fraction of gas than NW and RVWb. 
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Figure 9. Radial gas density profiles of galaxies found by FOF group finder at z = 1.98. It shows the gas overdensity as a function of 
comoving radius (or distance from the position of maximum density considered as the galaxy center), for four total (DM + gas + star) 
halo mass (Mhaio/A^©) ranges: lO'' - 10^° (top-left), lO^'' - lO" (top-right), lO" - lO^^ (bottom-left), and 10^^ - lO" (bottom-right). 
All the halos within each mass range (the number mentioned in the panels) are stacked over for each run, and the plotted curve denotes 
the median value in a radial bin. The grey shaded area enclose the 70th percentiles above and below the median in run RVWa (red 
curve), showing the typical radial scatter. The vertical dashed line is the analytical expression for the halo virial radius -R200 in comoving 
coordinates (from Eq.O, for Mhajo/A'^Q = 3 X 10^, 3 X 10^", 3 X 10^^, and 3 X 10^^ in the panels from top-left respectively. Details are 
given in t|3.6.1l 
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Figure 10. Gas temperature radial profile of galaxies at z = 1.98, in a format similar to Fig. [9] Effective temperature is shown for the 
multiphase gas particles ( i|3.6.2l l. 
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Figure 11. The gas carbon metallicity radial profile of galaxies at 2 = 1.98, in a format similar to Fig. (9] It shows the ratio of carbon 
mass fraction in gas to that of the Sun {Zc q = 0.002177), as a function of comoving radius, for four halo mass ranges. The dashed 
curves represent median- for the enriched particles only in each radial bin, i.e. not counting particles with zero-Zo ( ^3.6.311 . 



The dashed median Zc in Fig. [TT] is indistinguishable 
from the solid median Zc in the inner parts, since most 
of the gas is enriched. While at large r the dashed median 
Zc is higher, because the contribution of the non-enriched 
{Zc ~ 0) particles reduces the solid median Zc, as can be 
seen in the first three panels. Especially in the top-right 
panel the enriched median Zc's are 100 — 5000 times higher 
than the total median Zc's. 

Some features in Fig. [TT] are similar to the gas den- 
sity profiles f i|3.6.1|) . because metals are produced during 
SF which occurs in dense regions. At r < (0.4 — 0.6)i?2oo all 
the profiles show decreasing Zc going outward from center, 
with varying r-dependent negative slopes, which we later fit 
by a second order polynomial. The Zc profiles start to rise 
again from r ^ (0.4 — 0.6)-R200 and fall at larger r, form- 
ing a local peak at r > R2oo- Such a trend is visible in 
the top two panels for both solid and dashed curves, and 
in the bottom-left panel for the dashed curves only. It oc- 
curs because of a combination of reasons: the presence of 
surrounding substructures where more metals are produced 
in-situ by ongoing SF, and the spreading of metals by wind 
from the central SF regions into the CGM. These simulation 
results are consistent with observations which show breaks 
(changes of slope) in the radial metallicity profiles, and/or 
risin g metallicity gradients in the outer regions of galaxies 
(e.g. , IScarano fc Lepin3l2012l ) . 

In the inner r < (5 - 7)h~^ kpc CW and RVWa pro- 
files have a lower Zc than NW and RVWb, because wind 
suppresses central SF and transports some metal out. The 
trend reverses in the outer r > (5 - 7)h-^ kpc: CW, RVWa 
and RVWb produce a metallicity about 20 — 30 times higher 
than NW, because of accumulation of metal-enriched gas ex- 



pelled by wind. The differences are most prominent in the 
lower halo mass ranges (10^ — 10^" and 10^" — 10^^, top 
two panels), and decreases at higher masses. Massive halos 
of 10^^ — lO'^^ (bottom-left) still show significantly different 
profiles at large-r, while the profiles of the four wind models 
in the most-massive halos (bottom-right) are very similar, a 
reason for which is only the region inside -R200 is plotted for 
them. 

Model RVWb produces noteworthy changes in the Zc 
profile compared to NW: metallicity in the central part is 
slightly lower in RVWb reaching about half of NW, but 
starts to becomes larger at r > 5h^^ kpc up to 10 times 
higher than NW. This shows that the wind in RVWb, 
though least effective in other aspects, is substantially ef- 
fective in transporting metals away from SF regions into 
lower-density surroundings. 

We perform polynomial fits to the Zc radial profiles of 
runs NW and RVWa, within certain radial limits, for the 4 
halo mass ranges plotted in Fig. 1111 The median Zc (solid 
curves) versus r data are fitted with a second order polyno- 
mial of the form, log{Zc/Zc,e) = A + B(logr) + C(logr)^ 
and the resulting fit coefficients are listed in Table [S] The 
results for Mhaio = (10^° - 1O")M0 fitted between r = 
(1 - 40)/i"^ kpc are: 

log ( I = 0.0325 + 1.74(log r) - 3.10(log rf NW, 

V Zc,Q J 

= -0.194 + 0.542(log r) - 1.58(log rf RVWa. 

We infer from Fig. [TT] that the CGM gas at galac- 
tocentric distances close to and beyond -R200, within r ~ 
(30 - 300) ft" ^ kpc comoving, around galaxies of masses 
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A^haio/Af© — 10^ — 10^^, can give the best Zc observational 
diagnostic to distinguish between different galactic outflow 
models. 



3.64 CIV Fraction 

We compute the fraction of triply ionized carbon, CIV, of 
the gas particles in post-processing using phot oionization ta- 
bles d erived from CLOUDY (last described bv lFerland et alj 
1 19981 ). The relevant quantities used are redshift, gas density 
and temperature from simulation s napshot, and an ionizin g 
background from the HMOS tables (|Haardt fc Madaull200ll ) . 
as included in version 07.02.00 of CLOUDY. 

Fig. [T2] shows the radial profiles of the gas CIV frac- 
tion, /civ. There are varying positive slopes at small-r 
where /civ increases with radius, /civ attains a peak at an 
intermediate-r, which occurs at (20 — 60)/i~^ kpc depending 
on halo mass and wind model. Further out /civ decreases 
from its peak at larger-r, and have a negative slope up to 
300/i"^ kpc. 

In the top two panels, runs CW and RVWa produce 
higher /civ than NW and RVWb; the differences are signif- 
icant at r 7h~^ kpc reaching lO'^ — 10^ times, very small 
at larger r up to the peak of /civ, and increases again to 
10 — 70 times at r > 2Qh~^ kpc. The small-r positive-sloped 
regions in the bottom two panels show small differences: in 
the bottom-left RVWa has higher /civ than NW and RVWb, 
which are in turn higher than CW; whereas in the bottom- 
right NW has higher /civ than CW, RVWa and RVWb. In 
the bottom-left panel /civ in runs CW and RVWa reaches 
a ~ 50 times higher peak at a larger-r than runs NW and 
RVWb. The large-r negative-sloped regions shows the same 
trends as the top 2 panels, runs CW and RVWa producing 
higher /civ than NW and RVWb. 

We perform polynomial fits to the /civ median radial 
profiles, in a similar way as in i]3.6.3l but with 2 first-order 
polynomials of the form: log /civ = A + _B(logr), within 2 
radial limits. The resulting fit coefficients are listed in Ta- 
ble [S] The polynomials for Mhaio = (10^° - 10")Mq, fitted 
between r = (1 — Q)h~^ kpc are: 



log /civ = -7.01 + 6.01(logr) NW, 

= -3.10 + l.OO(logr) RVWa; 

and fitted within r = (6 — 2Q)h^^ kpc are: 

log /civ = -5.75 + 3.72(logr) NW, 

= -5.18 + 3.44(logr) RVWa. 



(12) 



(13) 



The overall behavior of /civ arises from the combined 
effects of density and temperature radial dependence of the 
gas (from Figs. and llOp in the ionization tables, with the 
major role played by the temperature. 

As a prediction for CIV fraction observations, we infer 
from Fig. [T5] that the inner gas at galactocentric distances 
r < (4 — 5) ft" ^ kpc comoving, around galaxies of masses 
Mhalo/M0 = 10^° - 10", can most-effectively distinguish 
between strong-wind and no-wind cases. 

3.7 Metallicity-Density Relation 

The carbon metallicity as a function of gas overdensity is 
plotted in Fig. 1131 a,t z = 1.98 in the top row, for the LB 



runs labeled by the color and plotting symbol. The larger 
boxsize is selected again for increased statistics, since more 
number of halos and more massive halos form in it. Gas par- 
ticles in three temperature bins are shown: 10^ — 10^ K in 
the top-left panel, 10^ — lO'^ K in the top-middle, and all 
the gas in the top-right. The solid curves denote the median 
Zc in (5-bins for each run. The grey shaded area enclose the 
70th percentiles above and below the median in run RVWa 
(red curve), showing the scatter at a given density. Follow- 
ing a format similar to Fig. Illl the solid curve medians and 
percentiles are computed considering all (both enriched and 
non-enriched) gas particles, while the dashed curves rep- 
resent median- for the enriched particles only i.e. those 
having Zc > 0. Note that if ^ 50% of the particles are not 
enriched, then the solid median would be zero. 

The dashed median Zc in Fig. 1131 top row, coincide 



with the solid median Zc at high densities, 5 ^ 10 



10* 



because most of the dense gas is star forming, generating 
metals and hence enriched. While at 5 < 10"^ the dashed 
median Zc's are several orders of magnitude higher than the 
solid ones (most clear in the top-left and top-right panels), 
implying that most or all of the low-density and underdense 
gas is not enriched. 

As an exception, the underdense gas having 0.1 ^ 5 ^ 2 
in the warm-hot 10^ — 10^ K phase (top-middle panel) is sig- 
nificantly enriched by the winds in runs CW and RVWa, 
so that all the particles (solid curves) have a median of 
Zc/Zc.Q ~ 0.05 — 0.08, just few times below the enriched- 
only dashed medians. Further trends discussed below ana- 
lyzes the dashed medians wherever the solid medians are 
below the plotting range. 

Overall the metallicity-density relation shows a nega- 
tive correlation at small-5 and positive correlati on at large-J. 
Simila r [Zc — S] relations were also found bv Tornatore et ahl 
(|2010t ) for the average metallicity of the warm-hot IGM at 
z — 0. We find that the threshold density for the slope 
turnover depends on temperature bin and wind model, lying 
between 5 = 30 — 1000. As a special feature, run CW has a 
slight increase in Zc from the smallest-!? to 5 ~ 10 in all the 
temperature bins, while run RVWa first decreases and then 
remains constant over the same range of densities. 

The wind models produce significantly different Zc in 
the low-density IGM at 5 < 10^ - 10^, with runs CW and 
RVWa enriching up to few 100 times higher than NW and 
RVWb. The metallicity induced by the four runs becomes 
similar at 5 > 10* in regions composed of SF gas. The warm 
phase of 10'' — 10^ K (top-left) presents the largest differences 
between the wind models: CW enriches the 5 ~ 30 gas about 
10 times more than RVWa, which is 10 times more than 
RVWb, which in turn is 50 times more than NW. 

RVWa produces higher Zc than CW at 5 < 5 (top-left 
panel), but for larger-5 values CW enriches more. The same 
trend (RVWa enriching more than CW at small S) is visible 
in the other two panels (top- middle and top-right), but the 
resulting metallicity becomes almost comparable in RVWa 
and CW at larger densities. 

Ap plying the pix e l opti cal depth technique to quasar 
spectra, ISchave et al.l l|2003l ) found a positive gradient of 
observed CIV metallicity with overdensity, measuring a me- 
dian [C/H] oc 0.65 log J, in the range log 5 = [—0.5, 1.8] and 
z = 1.8 — 4.1. Comparing with our Fig. 1131 dashed curves in 
the top row, some of the wind runs show a weak similar trend 
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Figure 12. Gas CIV (triply ionized carbon) fraction radial profile of galaxies at 2 = 1.98, in a format similar to Fig. [9] It shows the 
fraction of CIV (computed from simulation redshift, gas density, temperature and assuming a photoionizing background using CLOUDY 
ionization tables), as a function of comoving radius. More details in i|3.6.4l 




Pgas / < Pb > Pgas / < Pb > Pgas / < Po > 

Figure 13. Carbon metallicity versus overdensity of gas at 2 = 1.98 (top row, described in i|3.7|l . and at an earlier epoch z = 3.98 
(bottom row, discussed in i|3.8ll . for the LB runs labeled by the color and plotting symbol. Gas in three temperature ranges are shown: 
10" - 10^ K (left panels), 10^ - 10^ K (middle), and all the gas (right). The plotted solid curves denote the median value in each 
overdensity bin. The grey shaded area encloses the 70th percentiles above and below the median in run RVWa (red curve), showing the 
typical scatter. These medians and percentiles are computed using all (both enriched and non-enriched) gas particles. The dashed curves 
show the median- for enriched {Zq > 0) p articles only in overd ensity bins. The black dotted line in each panel shows the <5-range and 
\Zq — S\ slope obtained from observations by ISchave et al. I hoot ). 
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in this 5 range. In all the top panels of run CW and top- 
middle panel of RVWa, Zc show a weak positive correlation 
with 5, but with a slope ~ 4 times smaller than observed. On 
the other hand, runs NW and RVWb, top-left and top-right 
panels of RVWa, along with the all-particle solid medians of 
runs CW and RVWa in the top-middle panel, show a nega- 
tive correlation of Zc with 5, or almost constant Zc- At the 
same tim e, this negative correlation is con sistent with the 
results of iBaraiT Martel fc GermainI l|201lf) . Studying IGM 
enrichment from anisotropic AGN outflows, they found that 
at 2 ^ 2 the underdense regions (5 < 1) are enriched to 
higher metallicities, and the resulting [O / H] decreases with 
increasing IGM density, a trend more prominent with in- 
creasing anisotropy of the outflows. 

We perform polynomial fits to the median Zc (solid 
curves) versus 5 data of runs NW and RVWa, in a similar 
way as in i]3.6.3l and i]3.6.4l fitting with 2 first-order poly- 
nomials of the form: \og{Zc / Zc,q) = A + B{\og5), within 
2 overdensity limits. The resulting fit coefficients are listed 
in Table [4] for the 3 temperature ranges plotted in Fig. 1131 



The polynomials in the warm- hot (10^ 
within 10^ s; 5 < 10* are: 



10') K gas, fitted 



log 



/ Zc 



-3.81 + 0.656(log5) NW, 



= -3.22 -I- 0.546(log (5) RVWa; 
and fitted between 10* ^ 5 < 10^ are: 
Zc 



(14) 



log 



Zc,{ 



= -3.55 0.643(log 5) NW, 
= -2.70-F0.456(log(5) RVWa. 



(15) 



3.8 Redshift Evolution: Metallicity at z ~ 4 

We investigate the redshift evolution of two metallicity cor- 
relations: Zc versus r, and Zc versus 5, of the LB runs. 

The gas carbon abundance radial profiles of galaxies at 
an earlier epoch z — 3.98 are plotted in Fig. 1141 left panel, 
for halos in the mass range Mhaio/M© = 10^° - 10", having 
number of halos between A'haio = 4921 — 5336 for the four 
wind models. It is compared to the Zc profiles of the same 
halo mass range at z — 1.98 in the right panel of Fig. 1141 

We see that several of the median Zc versus r correla- 
tions are similar between z — 3.98 and z — 1.98, only the 
absolute metallicity values are lower at the earlier epoch. 
All the runs at 2: = 3.98 have Zc decreasing with r at 
r ^ (0.5 — 0.7)i?2oo, with r-dependent negative slopes. The 
trend of Zc profiles rising again at larger r values (because 
of in-situ metal generation in surrounding structures and 
spreading of metals by wind from the central SF regions) is 
weaker a.t z — 3.98; it is only visible for the dashed curves 
in the left panel. 

The differences between the wind models are smaller 
at 2: = 3.98 than z = 1.98, however there are signatures 
of suppression of central SF by the impact of winds, which 
transport metals out and accumulate them in the CGM. 
Analyzing the solid curves of the left panels, runs CW and 
RVWa have few times lower Zc than NW and RVWb in 
the inner r < (5 — 7)h^^ kpc. The trend reverses in the 
outer r > (5 - 7)/i"^ kpc; CW, RVWa and RVWb have few 
times higher metallicity than NW. The enriched-only dashed 



curves show a feature different from their solid counterparts 
for large r values, where they detach from the solid curves 
and increases with r; Zc in runs NW and RVWb are ~ 2 — 3 
times higher than CW and RVWa. 

The bottom row of Fig. [13] shows the gas carbon metal- 
licity versus overdensity at 2: = 3.98, for 3 temperature 
ranges, following the top row. Analogous to [Zc — r], sev- 
eral of the median Zc versus 5 correlations at z — 3.98 are 
similar to those a.t z — 1.98, with the absolute metallicity 
values being lower at earlier times. The differences between 
the wind models are smaller at z = 3.98 than z = 1.98. 

The underdense gas with 0.1 ^ 5 ^ 0.8 in the warm- 
hot 10^ — lO'^ K phase (bottom-middle panel) is already 
enriched by the winds in run CW, making the all-particle 
(solid curve) median Zc/Zc,q ~ 0.002 at z = 3.98, which is 
40 times lower than at z — 1.98. However, the enrichment 
at similar 5 values caused by the winds of model RVWa at 
z — 3.98 is substantially lower than that of CW by 1000 
times or more. 

Further trends at z — 3.98 discussed below analyze the 
dashed medians wherever the solid medians are below the 
plotting range. Starting from the lowest-(5, the metallicity- 
density in all the outfiow runs (except RVWb in the bottom- 
middle panel) show a shallow positive correlation at 5 ^ 30. 
At 2 = 1.98, such a feature is only present in run CW. 
This positive-sloped [Zc — S] relation that we obtain at z = 
3.98 in the w i nd m odels is consistent with observations by 
ISchave et al.l l|2003t ) who found a positive gradient of CIV 
metallicity at similar overdensities. 



4 SUMMARY AND CONCLUSION 

We explore new models of galactic winds performing 
hydrodynamic simulations using the TreePM-SPH code 
GADGET-3, and analyze their impact on the properties of 
the CGM at 2 = 2 — 4. Our outfiow implementation imparts 
kinetic feedback, in the energy-driven formalism, where the 
wind velocity has a positive correlatio n with galactocentri c 
radius v^ir) as seen in observations bv lSteidel et al.l l|201Cf ). 
We further investigate a halo mass dependent parametriza- 
tion of t he ra dially- varying wind, following observations by 
iMartinI l|2005l ). The simulations include additional subgrid 
physics: metal-dependent radiative cooling and heating in 
the presence of photoionizing background radiation; star for- 
mation; stellar evolution and self-consistent chemical enrich- 
ment using a fixed stellar IMF. 

The implementation of the new wind models in the code 
involves finding the distance of gas particles from their host 
galaxy center. We identify galaxies by running a FOF group 
finder on-the-fiy within a simulation at intervals of 1.001 
times the scale factor, and find stellar groups of at least 
32 particles, by linking over stars as the primary particle 
type, using a linking length 3 times smaller than that for 
obtaining DM halos. The location of the member gas particle 
with maximum density is considered as the galaxy center. 
Multiphase gas particles undergoing SF are stochastically 
selected and kicked into wind by giving their speed a one- 
time fu, boost in a direction perpendicular to the galaxy 
disk. Wind particles are also temporarily decoupled from 
hydrodynamical interactions. 

We simulate two different cosmological volumes (in or- 
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Figure 14. Gas carbon mctallicity radial profile of galaxies at an earlier epoch z = 3.98 in the left panel, compared to that at z = 1.98 
in the right panel, for halos with masses My^^i^/ Mq = 10^" — 10^^. Each panel shows the ratio of carbon mass fraction in gas to that 
of the Sun, as a function of comoving radius, in a format similar to Fig. 1111 but for one halo mass range. The dashed curves represent 
median Z^; values for the enriched particles only in each radial bin. See N3.8I 



der to increase the statistics) a smaller box of (5ft ^ Mpc)"^ 
comovme with 2 X 128^ DM and gas particles, and a larger 
box of (25/i~^ Mpc)^ with 2 x 320^ particles; using the fiat 
ACDM concordance model. For each volume, we perform 4 
runs investigating different galactic wind models: 
NW: no wind; 

CW: energy-driven wind with constant = 400 km/s; 
RVWa: radially varying wind with fixed parameters; 
RVWb: RVW with parameters dependent on halo mass. 

The main analyses of our simulations refiect the follow- 
ing processes: the outfiows expel gas away from the star- 
forming galaxies and suppress the SF; the gas is also metal- 
enriched; thus winds carry metals out and accumulate them 
in the CGM and IGM, enriching these lower density regions 
with metals. Our results are summarized below. 

- Outflow speed. 

The outflow gas velocity magnitude as a function of galac- 
tocentric distance of multiple wind-phase gas particles, ob- 
tained in our simulations, follows the given input subgrid 
wind speed, constant or varying with radius, in agreement 
with observations. 

- Star formation rate. 

Galactic wind feedback quenches SF at z < 8. In the 
(25/1^^ Mpc)^ runs at z ~ 2, the global SFRD is 4 times 
smaller in RVWa and 2 times smaller in CW than NW. 
At 2 ^ 5 RVWa causes a greater suppression than CW, 
and produces 2 — 4 times lower SFRD than CW aX z = 2. 
The SFRD at z ~ 4.5 - 10 is systematically (up to 2 - 10 
times) higher in the simulations, compared to observations. 
At later epochs z ^ 2 — 4.5, most of the observations (e.g., 
ICucciati et"alll2012l ) lie within the ranges of SFRD produced 
by the different wind models. 

- Gas and stellar mass functions and mass fractions. 
The gas mass function of galaxies in the {25h~^ Mpc)"^ 

runs aX z — 2.23 has the same slope between Mgas ~ (10^ — 
1O")M0, but in CW and RVWa is shifted leftward with 
respect to NW and RVWb by (2 — i)MQ\ because winds 
expel gas, causing a smaUer number of objects with high 



gas masses. Runs CW and RVWa produce a steeper stellar 
mass function than NW and RVWb and all of them are 
steeper than that in semi-analytical models and observations 
(e.g.. iBower. Benson fc Grain 2012'). There is a small excess 
of galaxies in RVWa compared to CW between Mstoiiar ~ 
(2 X 10* - 1O^)M0. 

In run NW the gas mass fraction decreases monotonically 
with halo mass; outflows (runs CW, RVWa, RVWb) flatten 
the Mgas /Af halo trend making it oscillatory. The stellar mass 
fraction in galaxies is the largest in run NW, and increases 
with Afhaio, such that NW and RVWb have the steepest 
increase, followed by CW, and RVWa is the flattest. 

- Thermal state of the gas. 

The [T — p] phase diagram shows that the outflow, soon 
after leaving the dense SF phase, goes through a cold (< 10* 
K) phase, due to the hydrodynamic decoupling. The thermal 
properties of the IGM look qualitatively similar in the 4 
runs. Model RVWa has a higher fraction of underdense to 
slightly overdense {S ~ 0.1 - 10), warm-hot (T ~ lO"* - 10*^ 
K) gas, than CW; both of them are significantly higher than 
NW and RVWb. 

- Gas kinematics. 

Projection of gas kinematics around the center of a galaxy 
of total halo mass Afhaio ^ 2 x lO^^M© at z = 2.12, shows 
that for NW and RVWb the gas outflows along the direction 
of least resistance of a low-density void. Run CW has most 
of the outflow inside r ^ 50 kpc forming a metal-enriched 
gas sphere, = 400 km/s not being enough to drive the 
wind to larger distances. The Uu,(r) in run RVWa produces 
an extended diffuse enriched gas outflow propagating per- 
pendicular to the galaxy disk (as expected from the model 
input), escaping to r > 100 kpc. The formation of a galaxy 
disk is visible in runs NW, RVWa and RVWb, where inside 
the {200h~^ kpc)^ projection volume about half of the gas 
is inflowing and the other half outflowing. Whereas the CW 
galaxy looks irregularly shaped, and more gas (64%) is un- 
dergoing infall than outflow (36%). The galaxy disk in run 
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RVWa is the largest in size and further investigations are 
needed to check if (r) can produce realistic disk galaxies. 

- Radial profiles. 

The radial profiles of gas properties around galaxy centers 
at 2 = 1.98 show most prominent differences between the 
models for lower halo masses (Mhaio/Mo = 10^ - 10^° and 
10^'' — 10^^), and almost uniform results at higher masses. 
The density and carbon metallicity profiles often form a lo- 
cal peak at r > -R200 of galaxies, because of the presence 
of smaller halos and surrounding substructures where more 
metals are produced by in-situ SF, and the spreading of 
metals by wind. 

• Gas density profile. Runs CW and RVWa have a 
lower overdensity in the inner r ^ lO/i"^ kpc, by 10 - 30 
times, than the NW and RVWb cases. While in the outer 
r > R200 , CW, RVWa and RVWb show a higher density 
than NW. 

• Temperature profile. The wind runs CW, RVWa and 
RVWb present a bimodal temperature distribution at 
r ^ lOh^^ kpc, composed of hot multiphase star form- 
ing gas and cold winds. Between (10 — 300)h~^ kpc, the 
NW model cause the highest-T (~ 10^ - lO'' K), followed 
by RVWb, then RVWa and CW. 

• Carbon metallicity profile. Runs CW and RVWa have 
a lower Zc than NW and RVWb in the inner r < 
(5 — 7)h~^ kpc, while in the outer parts runs CW, RVWa 
and RVWb produces higher Zc by 20 — 30 times than 
NW. Metallicity in RVWb becomes larger than in NW at 
r > 5h~^ kpc and up to 10 times higher at r > -R200- 
We perform second order polynomial fits to Zc versus r, 
whose coefficients are listed in Table (2] 

The carbon enriched {Zc > 0) gas fraction is 1 in the 
inner regions of galaxies; it reduces to (0.5 — 0.8) at -R200 
and to (0.45 - 0.7) at 300 h'^ kpc. Runs CW and RVWa 
enrich a higher (up to 1.4 times) fraction of gas than NW 
and RVWb. 

• CIV fraction around galaxies. CIV fraction profiles 
have varying positive and negative slopes at different r, 
attaining a peak between (20 — 60)/i~^ kpc. Runs CW and 
RVWa produces higher /civ than NW and RVWb. We 
fit the /civ radial profiles with 2 first-order polynomials, 
and list the coefficients in Table [31 The results for Mhaio = 
(10^" - 1O")M0, fitted between r = (1 - 6)h~^ kpc show 
the largest difi'erence: /civ cx: for NW, and /civ oc r^ 
for RVWa. 

• Observational Predictions. Inferred from our simula- 
tions, we predict that Zc observations of the CGM gas at 
galactocentric distances in the range r ~ (30 — 300) 
kpc comoving, around galaxies of Afhaio = (10^ — lO^^)Af0, 
can best distinguish between different galactic outfiow sce- 
narios. And CIV fraction observations of the inner gas in 
the range r < (4 — 5)ft~^ kpc comoving, around galaxies of 
Aihaio ~ (10^" — lO^^)Af0, can most-effectively distinguish 
between strong-wind and no-wind cases. 

- Metallicity as a function of density. 

The underdense 5 = 0.1 — 2 warm-hot 10''' — 10^ K phase is 
significantly enriched by the winds in runs CW and RVWa, 
so that the all-particles median is Zc/Zc,q ~ 0.02—0.08. Zc 
versus 6 at z = 1.98 shows a negative correlation at small-5, 
and positive correlation at large-5. The only exception is run 
CW where Zc has a small increase between (5 ~ 0.1 — 10, 



while Zc in RVWa first decreases and then remains fiat. The 
low-density IGM with 5 < 10^ — 10'^ is significantly enriched 
in runs CW and RVWa, up to few 100 times more than NW 
and RVWb. In the warm 10* - 10^ K phase, CW enriches 
the 5 ~ 30 IGM 10 times more than RVWa, which enriches 
10 times more than RVWb, which in turn enriches 50 times 
more than NW. 

We perform 2 first-order polynomial fits to the [Zc — S] 
correlation, whose coefficients are listed in Table 3] As an 
example, for the warm-hot (10^ — 10^) K gas, the fits within 
10^ ^ 5 < 10'' are: Zc oc 6° '^^'^ for NW, and Zc oc 5°'^'"' for 
RVWa. 

- Redshift evolution. 

Several of the [Zc — r] and [Zc — 5] correlations at 2: = 
3.98 are similar to those at 2 = 1.98, with lower metallicity 
values and smaller differences between the wind models at 
the earlier epoch. The trend of Zc radial profiles rising again 
at large-r is almost absent at z = 3.98. For Mhaio/A/© = 
IqS _ 10^0 and 10^° - lO", the enriched-only median Zc in 
runs NW and RVWb at large-r are 2 — 3 times higher than 
in CW and RVWa. 

The 2 = 3.98 underdense 5 = 0.1 - 0.8 warm-hot lO'' - lO'' 
K gas is enriched by the winds in run CW (but not in RVWa) 
making the all-particle median Zc/Zc,q ~ 0.002, which is 
40 times lower than at 2 = 1.98. The [Zc — 5] relation at 
2 = 3.98 in all the wind runs (except RVWb in the warm-hot 
phase) shows a shallow positi ve correlation over 5 ~ 0.1 — 30, 
consistent with observations jSchave et al.|[20o3 ). 

In summary, we have found that the wind model 
with the radially varying velocity dependent on halo mass 
(RVWb) is the least effective in modifying IGM related prop- 
erties, with results similar to the no- wind (NW) case, except 
that it substantially enriches the low-density CGM. 

The impact of the model RVWa, for which the veloc- 
ity is increasing as a function of galactocentric distance, is 
instead similar to the energy-driven constant-velocity im- 
plementation CW. However, it shows interesting differences 
that deserve to be further investigated. RVWa causes a 
greater suppression of SFR than CW at 2 ^ 5, this could 
have implications for the galaxy downsizing scenario. RVWa 
also produces galactic disks larger than all the other wind 
models, and one can study if the radially varying outflow for- 
malism can generate more realistic disk galaxies. Run RVWa 
has a higher gas fraction than run CW in the low-density 
{5 ~ 0.1 - 10) warm-hot (10'' - 10'' K) phase of the IGM, 
which could shed light on the missing baryon problem. 

We see different trends of Zc versus 5 between the CW 
and RVWa models at 5 ^ 10, and CW outflows generally 
produce a higher and earlier enrichment of some IGM phases 
than RVWa. To explore such IGM metal-enrichment differ- 
ences, future progress in this field should include computing 
more observable statistics from the simulations, e.g., Lyman- 
a fiux, simulated quasar spectra, and compare them with 
observations of CGM and IGM at different impact parame- 
ters from galaxies. 
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Table 1. Simulations Parameters ^ 



Run Lbox ^ NpB,rt T^gas Lsoit ° Galactic Wind Feedback 

Name [ft"^ Mpc] [h~'^Me] [/i"^ kpc] 



Smaller-Box Runs : SB * 
NWt 5 2 X 128^ 7.66 x 10^ 0.98 No Wind 

CWt 5 2 X 128^ 7.66 x 10^ 0.98 Energy-driven constant-velocity v = 400 km/s 

RVWat 5 2 X 128^ 7.66 x 10^ 0.98 Radially varying with fixed parameters ^ 

RVWbt 5 2 X 128^ 7.66 x 10^ 0.98 RVW with halo mass dependent parameters ^ 



Larger-Box Runs : LB ' 
NW 25 2 X 320^ 6.13 x lO'* 1.95 No Wind 

CW 25 2 X 320^ 6.13 x 10® 1.95 Energy-driven constant-velocity v = 400 km/s 

RVWa 25 2 x 320^ 6.13x10*^ 1.95 Radially varying with fixed parameters 

RVWb 25 2 X 320^ 6.13 x lO'' 1.95 RVW with halo mass dependent parameters ^ 



^AU simulations have the same physics described in §2, with the wind model varied. 
^ Lhox ~ Comoving side of cubic simulation volume. 

A'^part = Number of gas and DM particles in the initial condition. 

rrigas ~ Mass of gas particle (which has not undergone any star- formation). 

Lsoit ~ Gravitational softening length (of all particle types). The minimum gas smoothing length is set to a 
fraction 0.001 of Z/soft- 

'Run names ending with "t" are smaller boxsize runs SB, stopped at z ~ 2. 

^Parameters of radially varying wind model ( ^2.3^ : rmin ~ lh~^ kpc, Rch ~ WOh~^ kpc, Umax ~ 800 km/s, 
Q = 1.15. 

'^Parameters of radially varying wind model dependent on halo mass ( i]2.3p : t^max ~ 2^/ GAfhaio/-R200, and 

RaB ~ R200- 

'Runs with larger boxsize LB, stopped at z ~ 2. 



Table 2. Polynomial Fit of Carbon Metallicity Radial Profile 



-^^halo,min -^^halo,max ^min 


' max 


Run 


Fit Coefficients 


[Mq] [Mq] [h-^ kpc] 


[h-'^ kpc] 


Name 


ABC 



1Q-' 10'" 1 30 NW 0.0580 0.667 -2.37 

10^ 10^° 1 30 RVWa -0.427 0.891 -1.77 

10'° 10" 1 40 NW 0.0325 1.74 -3.10 

lO"' lO" 1 40 RVWa -0.194 0.542 -1.58 

10" 10'^ 1 300 NW 0.338 0.478 -1.02 

10" 10'^ 1 300 RVWa -0.00471 0.640 -0.913 

10" 10'^ 1 300 NW 0.449 -0.147 -0.349 

10'^ 10'^ 1 300 RVWa 0.529 -0.309 -0.227 



^Coefficients from fitting median Zc versus r data with a second order polynomial of the form, 
log{Zc/Zc,Q) = yl-f -B(logr)-|-C(logr)^. Fits are calculated for 2 runs from Fig.[Tl]at z = 1.98, 
for the 4 halo mass ranges (4 panels in the figure) within radial ranges indicated here. 

''Minimum and maximum halo mass indicating the range of the median Zc result. 

"^Minimum and maximum radius from galaxy center within which the fit is done. 
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Table 3. Polynomial Fit of CIV Fraction Radial Profile 



A/f ^ 

-'^-^halOjinin 


-'^-^halojmax 


Run 


''^l,min 


^l.max 


Fit Coeffs-1 


r2,min " 


^2, max 


Fit Coefrs-2 


[Mq] 




Name 


[h-^ kpc] 


[ft-i kpc] 




Bi 


[h-^ kpc] 


[h-^ kpc] 






10^ 


10"' 


NW 


1 


6 


-3.94 


2.25 


6 


20 


-5.17 


3.60 


10^ 


10^" 


RVWa 


1 


6 


-2.69 


0.642 


6 


20 


-4.65 


3.22 


10^° 


10" 


NW 


1 


6 


-7.01 


6.01 


6 


20 


-5.75 


3.72 


10^° 


10" 


RVWa 


1 


6 


-3.10 


1.00 


6 


20 


-5.18 


3.44 


10" 


10^2 


NW 


2 


7 


-11.96 


9.14 


7 


20 


-5.08 


2.15 


10" 


10^2 


RVWa 


2 


7 


-10.55 


7.98 


7 


20 


-5.97 


2.93 


10^2 


10" 


NW 


4 


10 


-17.05 


14.00 


10 


30 


-4.35 


1.07 


10^2 


10" 


RVWa 


4 


10 


-14.23 


9.74 


10 


30 


-7.63 


3.45 



''Coefficients from fitting median /civ versus r data witli two first-order polynomials of the form: log /civ = Bi (log r) 
within [ri,min, n.max], and log/civ = A2 + B2(\ogr) within [r2,min, r2,max]. Fits are calculated for 2 runs from Fig. [121 at 
z = 1.98, for the 4 halo mass ranges (4 panels in the figure). 

'^Minimum and maximum halo mass indicating the range of the median /civ result. 

'^Minimum and maximum radius from galaxy center within which the fit is done. 



Table 4. Polynomial Fit of Carbon Metallicity-Overdensity Relation ^ 



rp b 

mill 


rp b 

max 


Run 


5l,miii 


Ol ,max 


Fit Coeffs-1 


52,min 


c r c 
02, max 


Fit Coeffs-2 


[° K] 


r K] 


Name 






^1 


Bi 






A2 


B2 


10* 


10^5 


NW 


10^ 


lO'^ 


-7.54 


1.43 










10* 


lO'^ 


RVWa 


10^ 


lO'^ 


-5.59 


0.937 










10^ 


10^ 


NW 


10^ 


10* 


-3.81 


0.656 


10* 


10^ 


-3.55 


0.643 


10^ 


10^ 


RVWa 


10^ 


10* 


-3.22 


0.546 


10* 


10^ 


-2.70 


0.456 


10' 


10* 


NW 


10^ 


10^ 


-7.62 


1.54 


10^ 


10* 


-1.20 


0.247 


10^ 


10* 


RVWa 


10^ 


10^ 


-5.98 


1.19 


10^ 


10* 


-1.74 


0.300 


''Coefficients from fittiuf 


5 median Zc 


versus & 


= Ps!^/{pb) dati 


a with two first-order polynomials of 


the form: log(.Zc/^c,Q) = 


^1 + Bi(lo! 


l5) within [(5i,min, 


^1 ,max] , 


and log(Zc/^c,0) 


= A2 + B2{logS) 


within [82, 


mill, ^2, max]. FitS 


are calculated for 2 runs from Fig.ll3lat 


z = l.i 


38, for the 3 j 


jiven temperature 


ranges (3 panels in 


the figure). 

















''Minimum and maximum temperature indicating the range of the median Zc result. 



'^Minimum and maximum overdensity within which the fit is done. 
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